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Abstract 

We classify models generating large local-type non-Gaussianity into some cat- 
egories by using some "consistency relations" among the non-linearity parameters 
/nl ' ^nl an d 5n£ ' wn ich characterize the size of bispectrum for the former and 
trispectrum for the later two. Then we discuss how one can discriminate models of 
large local-type non-Gaussianity with such relations. We first classify the models by 
using the ratio of T^ al /(Gf^ 31 /5) 2 , which is unity for "single-source" models and 
deviates from unity for "multi-source" ones. We can make a further classification 
of models in each category by utilizing the relation between /^L al an d <?NL a '- Our 
classification suggests that observations of trispectrum would be very helpful to dis- 
tinguish models of large non-Gaussianity and may reveal the generation mechanism 
of primordial fluctuations. 
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1 Introduction 



Non-Gaussianity of primordial fluctuations is now one of the important observables to 
probe the physics of the early Universe. Although the inflaton, which acquires quan- 
tum fluctuations during inflation and generates (almost) scale-invariant and adiabatic 
primordial fluctuations, has been considered to be responsible for the origin of density 
fluctuations in the Universe, this scenario might not be the one realized in the nature. 
Primordial fluctuations generated from the inflaton during slow-roll inflation are expected 
to be almost Gaussian, however, current observations such as WMAP indicate that the 
primordial fluctuations might deviate from Gaussian ones and the deviation could be rel- 
atively large compared to the one expected from the "slow-rolling" inflation models. (For 
recent analyses, see [IHI].) The degrees of non-Gaussianity are often represented by a 
so-called non-linearity parameter /nl, which characterizes the size of bispectrum of the 
curvature perturbation. Depending on the momentum distribution of the bispectrum or 
the shape of three point function, three types of /nl have been discussed in the litera- 
ture [HE]: local, equilateral and orthogonal types. The limits on these /nls have been 
obtained as pQ: -10 < /^ al < 74 for the local type, -214 < f^ il < 266 for the equi- 
lateral type and —410 < f^l h ° e < 6 for the orthogonal type (95% C.L.). Although all 
of these are still consistent with Gaussian fluctuations, they may give some hint of non- 
Gaussian ones since the central value of /nl for some types is away from the zero. If the 
primordial fluctuations are confirmed to be deviated from Gaussian ones in the future, 
which may well be probed with more precise observations such as Planck [6], fluctuations 
from the "slow-rolling" inflaton would be excluded at least as a dominant mechanism of 
the generation of density fluctuations. However, many other mechanisms have also been 
discussed in the literature and quite a few of them can generate large non-Gaussianity. 

Now we have plenty of models generating large non-Gaussian primordial fluctuations. 
Thus the question we should ask next is "how can we differentiate these models?" In 
this paper, we discuss this issue by using bispectrum and trispectrum of the curvature 
perturbation. Although, as mentioned above, we can divide models into some categories 
depending on the shape of the three point functions (local, equilateral and orthogonal 
typesf^l the shape is not enough to differentiate models since there remain many models 
for each type of the shape. Furthermore, /nl predicted in those models can fall onto almost 
the same value by tuning some model parameters. Thus, obviously, the determination of 
/nl is not enough to pin down the models of large non-Gaussianity even if / N l is found 
to be large in the future. The purpose of this paper is to pursue the strategy of how 
one can differentiate models of large non-Gaussianity. To this end, we consider higher 



order statistics such as the trispectrum in addition to the bispectrumF 2 !. The size of the 



There are various inflationary models that can produce large /nl" > where scalar fields have some non- 
canonical kinetic terms [71413] . The orthogonal- type /nl' 105 approximately arises from a linear combination 
of higher-derivative scalar-field interaction terms [5j[12] . 
# 2 Another possible direction may be the one using the scale dependence of /nl [SHE] • 
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trispectrum can be parametrized by other non-linearity parameters tnl and on tF 3 ! and the 
importance of the trispectrum has been emphasized in some literature [2QM25] . However, 
here we give a systematic study of the bispectrum and the trispectrum of models with 
large non-Gaussianity and make some classifications by using the "consistency relations" 
between the non-linearity parameters /nl, tnl and #nl- 

Among three types mentioned above (local, equilateral and orthogonal types), we focus 
on local-type models in this paper. As will be shown, by using the "consistency relations," 
we may be able to discriminate models of large non-Gaussianity. 

The organization of this paper is as follows. In the next section, we summarize the 
formalism for the discussion of local-type non-Gaussianity of primordial fluctuations. Then 
in Section [31 we first make a classification of models by using some relation between the 
non-linearity parameters /nl ; tnl and #nl m a systematic manner. Some example models 
will also be discussed. The final section is devoted to summary and conclusion of this 
paper. 



2 Formalism 

First we summarize the formalism to discuss local-type non-Gaussianity of primordial 
fluctuations, which is expected to be generated from super-horizon dynamics of curvature 
fluctuations. For this purpose, we adopt the 5N formalism [26H29] . In this formalism, the 
super-horizon scale curvature perturbation on the uniform energy density hypersurface £ 
at some time t — tf is given by the perturbation in the number of e- folding measured from 
the initial time i* to the time tf. Here we take the initial time i* to be the time shortly 
after the horizon crossing and the initial hypersurface to be a flat one. Then the curvature 
perturbation is given, up to the third order, as 

C(t f ) ~ NJtf + ^N ab 6<p a M + ^NatJtfSrfdrf , (1) 

where a subscript a, b and c labels a scalar field, which is assumed to be Gaussian fluc- 
tuations 5ip a at t — t* in the following discussion, and N a = dN/dip^ and so on. The 
summation is implied for repeated indices. 

The power spectrum P^, bispectrum and trispectrum of the curvature perturba- 
tion are given by 

(^ 1 C A r 2 ) = (2vr) 3 P ? (fc 1 )5(fc 1 + fc 2 ), (2) 

(Ck^Cks) = WB c {h, k 2 , h)S(h + k 2 + k 3 ), (3) 

# 3 Current observational limits for r^£ al and 5 ^£ al are given by -0.6 x 10 4 < r^ 1 < 3.3 x 10 4 (95% 
C.L.) and -7.4 x 10 5 < g^£ al < 8.2 x 10 5 (95% C.L.) from cosmic microwave background observations [23] . 
The authors of gl] also obtained -3.5 x 10 5 < c^£ a! < 8.2 x 10 5 (95% C.L.) from large scale structure 
observations by assuming that f^£ al = and t^^ 1 = 0. 



2 



(C^C*^) = (2vr) 3 T ? (A; 1 , fa, fa, fa)5{fa + h + fa + fc 4 ), (4) 
where S( and can be written as 

B c (fa, fa, fa) = ^'° L cal (P^k^ifa) + P f (fc 2 )P c (fc 3 ) + P ( (h)Pd k i)) , (5) 

?c(A;i, fc 2 , fa, fa) = r£ L caI (P c (k 13 )P c (fa)P ( (fa) + 11 perms.) 

+ ^nT' (Pdh)Pdh)Pdh) + 3 perms.) , (6) 

with fa 3 = \ fa + A^3 j . Here f^ al ,T^ al and (?N£ al are non-linearity parameters of the local 
type. From Eq. ([1]), the power spectrum of the curvature perturbation is given by 

P c (k) = N a N a P s (k) , (7) 

where P$(k) is the power spectrum for fluctuations of a scalar field: 

(S^MkJ = i^) 3 ^ (h + h) Ps(h) = (2vr) 3 ^ (% + fa) ^-V s (fa), (8) 

with Vs(k) = [H^/2tt) 2 and being the Hubble parameter at t = t*. 

Since 5ip a are Gaussian fields, the leading order contributions to the bispectrum are 
given by four-point functions of 5ip a . The next leading order contributions are given by 
six-point functions. As is explained in detail in [30|[3Tj. we can systematically classify 
those contributions by assigning a diagram to each contribution. In terms of the diagram 
approach, the four-point function of 5ip a is represented by the tree diagram and the six- 
point function by the diagram containing a single loop. This can be easily generalized 
to higher order correlation functions of ( (and also to the power spectrum). The leading 
order contributions to the n-point function of ( are given by 2{n — l)-point function of 
5ip a , which are represented by the tree diagrams. The next leading order contributions are 
given by 2n-point function, which are represented by the one-loop diagrams. 

In many cases, the contributions from the loop diagrams are tiny and can be safely 
neglected. Therefore, we first drop the loop corrections and consider only tree contribu- 
tions. We will discuss a case in which effects from loop corrections cannot be neglected in 
a later section and an appendix. 

By calculating the bispectrum and the trispectrum at the tree level, we obtain the 
non-linearity parameters /nl^nl and #nl as 



6 .local _ N a N b N<* 

5 /NL (N c N<) 2 ' 1 J 



_ local N a N b N^N c b 
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54 local = N abc N*N b Nt 

25 5nl (N d N d ) 3 ' 1 ' 

In the following, since we concentrate on non-Gaussianity of the local-type models, we 
omit a superscript "local" unless some confusions arise. 



3 Models of local-type non-Gaussianity 
3.1 A classification of local type 

Even if we limit ourselves to models generating large local-type non-Gaussianity, there still 
remain many possibilities. To discuss how we discriminate those models, here we classify 
the local-type models into some categories. For this purpose, we first give a particular 
attention to the relation between /nl and tnl- Notice that, when a single scalar field a is 
only responsible for density fluctuations in the Universe, 5N formalism indicates: 

g/NL--^, 25 9 ^~^f ^) 

From these equations, we can easily find the following relation between /nl and tnl: 

6 ^ 2 



TNL = U/nlJ • (13) 

Since some simple models such as the pure curvaton and the pure modulated reheating 
scenarios fall onto this type, we use the above relation to make a classification. We 
call models in which Eq. (TT3~]) is satisfied as "single- source model." Notice that tnl is 
determined once /nl is given in models of this class. 

The next category is for models which do not have a universal relation between /nl 
and tnl, although there exists an inequality as given below. This is a general situa- 
tion where there are multiple sources of density fluctuations. Thus we call this category 
"multi-source model." An example of this type includes mixed fluctuation models, where 
fluctuations from both of the inflaton and another scalar field such as the curvaton can be 
responsible for density fluctuations. Since the inflaton gives almost Gaussian fluctuations, 
non-Gaussianity mostly originates from fluctuations of the other source in such a case. 

The third category would be for models where there is some universal relation between 
/nl an d Tnl- In fact, this definition includes models of the first class discussed above as 
a special case. However, the "single source model" is somewhat a very special type, and 
thus we treat them as a separate category. One of examples of this third category is so- 
called "ungaussiton" model [34TI36] . in which the relation tnl oc holds. In this model, 
the second order perturbation of the "ungaussiton" field dominates the linear order one. 
Since the dominance of the second order perturbations indicates that the fluctuations are 
completely non-Gaussian, we need another source which gives a Gaussian perturbation at 
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linear order. Thus this model would require multi-sources of fluctuations in this sense. 
Hence we call models in this third category "constrained multi-source model." 

Here it is worth noting that the following inequality holds in the /nl~tnl plane: 



which can be obtained by using Cauchy-Schwarz inequality and was first found by two of 
the present authors (T.S. and M.Y.) [H7]. The derivation of this inequality ( tHJ) in [H7] 
was based on the expressions for /nl and tnl given in Eqs. and f llOl) . which are valid 
at the tree level. In this paper, we further show that, even when some loop corrections 
are dominant in /nl and tnl, the inequality f fl4|) holds as far as loop contributions are 
subdominant in the power spectrum, which is required from current observations. Since 
all models generating large local-type non-Gaussianity known to date satisfy the above 
assumptions, the inequality f[T4"j) is very important to test the local-type models. Thus 
we call the inequality ( IT4|) "local-type inequality" in the following. We discuss this issue 
in detail again in Section 13.41 Also notice that single-source models correspond to the 
boundary of this local-type inequality. 

In the following sections, we discuss various models for each category. We start with 
"single-source model," then discuss "multi-source model" and "constrained multi-source 
model" afterwards. In Fig. [U some models of these categories are shown in the "Jni^Tnl" 
diagram, from which we can have some idea of to what extent models of each category 
can give different predictions on / NL and tnl- We also give a summary of the categories 
and their examples in Table [TJ 

Now we have categorized models of local-type into three classes by using the key quan- 
tity t nl /(6/nl/5) 2 . However, since each category still includes some (or many) possible 
models, we may need another quantity to discriminate them. For this purpose, we can fur- 
ther utilize the relation between /nl and #nl- Although the /nl^nl relation can change 
depending on the model parameters, we can roughly divide models into three types further 
by looking at their relative size. As we will argue in the following sections, the relation 
|<?nl| ~ |/nl| holds in some models, then we call such models "linear g NL type." In other 
models, <?nl could be suppressed compared to /nl, i-e. <7nl ~ (suppression factor) x / NIj , 
which we denote this type of models as "suppressed #nl type." The other type is "en- 
hanced #nl type" in which the relation between / NL and #nl can be given as #nl ~ /nl 
with n > 1 (but in most models discussed in this paper, n — 2). By using the /nl _ t"nl and 
/nl^nl relations, we may be able to discriminate models well. In Table [TJ the /nl^nl 
relations for example models are shown and depicted in Fig. [2j 

3.2 Single-source model 

If future observations confirm the relation r NL = (6/nl/5) 2 , models of large non-Gaussianity 
we should pursue would be the ones in this category. Models categorized in this class in- 




(14) 
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Figure 1: /nl^tnl diagram. The relation between /nl and tnl is shown for three categories: 
single-source, multi-source and constrained multi-source models. For multi-source and 
constrained multi-source models, the cases for some representative explicit models (mixed 
curvaton and inflaton, mixed modulated reheating and inflaton, and ungaussiton models) 
are plotted. All the three categories satisfy the inequality tnl > (6/nl/5) 2 as far as loop 
contributions are subdominant in the power spectrum. The region with tnl < (6/nl/5) 2 
is shaded with gray. 
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Category 


/nl-tnl relation 


Examples and /nl-^nl relation 


Single-source 


tnl = (6/nl/5) 2 


(pure) curvaton (w/o self-interaction) 
[SNL = -(10/3)/ NL -(575/108)] (a) 


(pure) curvaton (w/ self- interaction) 

T A £ 2 7~> X /~i 1 (^0 

[fi'NL = ^NQJnL + -DNQJNL + <-^NqJ 


(pure) modulated reheating 
[^nl = 10/nl-(50/3)] (c) 


modulated-curvaton scenario 

o 1/2 ,3/2] W 
fiNL - ^dcc/NL 


Inhomogeneous end of hybrid inflation 

[fiNL = (10/3)r/cr/NL] 


Inhomogeneous end of thermal inflation 
[ 5NL = -(10/3)/ NL -(50/27)] (e) 


Modulated trapping 
[g NL = (2/9)^ L ] (/) 


Mult i- source 


tnl > (6/nl/5) 2 


mixed curvaton and inflaton 
[gNL = -(10/3)(i?/(l + R))f NL - (575/108)(i?/(l + R)f} {9) 


mixed modulated and inflaton 
[fi'NL = 10(^/(1 + i2))/ NL - (50/3)(fl/(l + /?)) 3 ] W 


mixed modulated trapping and inflaton 
[fi'NL = (2/9)((l + R)/R)fL = (25/162)r NL ] (i) 


multi-curvaton 

[fiNL = C mc / NL , fi'NL = (4/15)/2 L ] (j) 


Multi-brid inflation (quadratic potential) 

[fiNL = -(10/3)T7/ NL , fiNL = 2/2 L ]( fc ) 


Multi-brid inflation (linear potential) 

[fi'NL = 2/nl] ' 


Constrained 
multi-source 




ungaussiton (C ~ 10 3 , n = 4/3) 



(a)p or j^c case w ;th r\j GC <C 1. 

^'AnQjBnq and Cnq are given in Eqs. ||23|I - H25( and this expression is for r^ cc <C 1. 
(^Facra = is assumed. 

( d ) This relation holds in the Region 2. For other cases, see text. 

(e) gin = o j s assumed. 

W'A = ct/M and m = ga are assumed. 

WA quadratic potential and r dcc -C 1 are assumed for the curvaton sector. R = PY 7 ' /P^ is the ratio of the power 

spectra. This relation can also be written as <?nl — — (24/5)(/^ l /tnl) — (9936/625)(/^ l /t^ l ). 

'^^Tctcto- = is assumed for the modulated reheating sector. This relation can also be written as 

9NL ~ (72/5)(/3 L /r NL ) - (31104/625)(/« L /r3 L ). 

(OA = a/M and m = go are assumed for the modulaton sector. 

( J )The former and the latter relations are for the cases where both curvatons are subdominant and dominant at their 
decay, respectively. C mc is 0(1) coefficient and always negative. 

( fe )The former and the latter relations are for the equal mass and the large mass ratio cases, respectively. 
(''For the equal mass case with gi = g2- 



Table 1: Summary of the categories and their examples. 
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Figure 2: /nl^nl diagram. The relation between g NL and /nl is plotted for models given 
in Table [TJ 
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elude the pure curvaton and the pure modulated reheating scenarios, and so on. In the 
following, we look at some of models in this class more closely. 



3.2.1 Pure curvaton model 

To make a prediction on the non-linearity parameters in the pure curvaton model [3"S1 - I40j#4 
first we need to specify the curvaton potential. Although in most study of the curvaton pTIT 
l5T] . a quadratic potential is adopted, it is also possible to have a self- interaction term. 
Thus we assume the following potential [2 H l22 |l52H55j ^l: 

n°) = Irnla 2 + \mt (j-J , (15) 

where m a is the mass of the curvaton, A is the coupling and n is the power for a self- 
interaction term. In the following discussion, we assume that the self-interaction term 
is subdominant compared to the quadratic one. In general, it is possible that a self- 
interaction term dominates first, and then the quadratic term becomes dominant at later 
time. However, in such a case, rigorous numerical calculations will be needed. Thus we 
restrict ourselves to the "nearly" quadratic potential case when we include a self-interaction 



term * 6 



In fact, to make a concrete prediction on non-Gaussianity in this model, we need to 
specify the decay rate T a , the mass m a and the initial value for the curvaton o~*. However, 
as far as the non-linearity parameters are concerned, the information on these quantities 
is encoded into the parameter r^cc, which corresponds to the ratio of the energy density of 



the curvaton to the total one at the decajej. The definition of rd GC and its approximate 

(16) 



expression using T a , m a and a* are given by 



'"dec — 



4prad + 3p CT 



dec 



M P Vr CT /m 



where M p \ ~ 2.4 x 10 18 GeV is the reduced Planck mass. Thus in the following, we do not 
specify explicit values for r CT , m a and cr*. 

In the curvaton model, adopting the sudden-decay approximation and assuming radiation- 



dominated background from the period of the start of its oscillation to its deca\r 8 I. the 



^ 4 We later also consider a mixed model where density fluctuations originate both from the curvaton 
and the inflaton. To be definite, we call the model where only the curvaton generates density fluctuations 
the pure curvaton model. 

# 5 For the analyses of other non-quadratic type potentials, see [56ll59j . 

# 6 For a general case, see |54[[55] . 

* 7 When the initial amplitude for the curvaton is large enough, the curvaton can drive the second 
inflation after the standard inflation. In this case, this statement does not hold true. However, the value 
of /nl is 0(1) in such a situation, hence we do not discuss such a case here. 

# 8 In Ref. [22], it has been shown that the background equation of state also affects the curvature 
perturbation generated from the curvaton, in particular, its non-Gaussianity. 
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curvature perturbation can be given bjcfl, 



i 2 ^"osc c ~^ 

Ccur 7T '"dec 00"* -\- — 
<J °~ osc 

4 

+ 81 



3r d « ( 1 + ) - 4rL - 2r dec 

9rdcc /' c7 osc lT osc 1 g ^osc^osc A _ g^2 A _|_ °° sc0 osc 



2 

osc \ / r _ \ 2 



0"n 



(50"*; 



CT /3 a f2 I dec I ~ /2 

^osc ^osc / \ '-'osc 



1 - 9^#^ ) + 10rL + 3rL 

^"osc 



/ \ 3 



0" n 



(oV*) 3 , (17) 



where <r osc denotes the value of the curvaton field at the time of the start of oscillation 
and the prime indicates a derivative with respect to cr*, namely a' osc = da osc /da^ and so 
on. o"osc and its derivatives represent the non-linear evolution of the curvaton after the 
horizon exit. 

From the expression of the curvature perturbation given above, the non-linearity pa- 
rameters /nl and #nl can be evaluated as 

5 /NL = 2^l, 1+ ^srJ- 2 - r "- (18) 

54 9 f(T^Cc ^o^wCA 9 / ffosc< sc 

^NL = T-2- 35 + 3 _,o - — 1 + 



25 4r 2 \ rr' 3 rr' 2 / V rr' 2 

*° ^' dec V u osc u osc / ' dec \ u os 

1 „eT n8c Cr" 



osc 



+ 1- 9^^ + 10r dec + 3rL- (19) 



When the potential has a pure quadratic form, we have 

^ = ^, <sc = 0, <C = 0. (20) 

Oosc 

Thus, in this case with r^cc < 1, we can find the following relation between /nl and 

10 575 

9NL = -y/NL ~ TJj^ + ( r dec)- (21) 

On the other hand, when we cannot neglect the contribution from a self-interaction term, 
the non-linear evolution of the curvaton field after the horizon exit changes the relation 
between /nl and <?nl, for r^ec < 1, as 

fi'NL — ^4nq/nl + -Bnq/nl + Cnq, (22) 



# 9 In the curvaton model, isocurvature fluctuations can also be generated. However, we do not consider 
such a case here. For the study of isocurvature fluctuations in the model, see [4Tll43ll49ll50j . 
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where the coefficients are given by 



9 / it 2 rr'" it it" \ / it it" 

Z ' °osc°osc i qC'osc'^osc I / 1 i °osc°osc 



A I OSC OSC i 9 Of- OSC 1 I 1 I ohL OSC | /no\ 

A m ~ o \ —Ts — + 6 ~ 7n — 1 + — 32 — J ' W 



S V rr' 3 rr' 2 I \ rr'' 

\ ^OSC OSC / \ "osc 

D 10 ft p-^osc^osc q^osc^oscA ( -I | ^oscCoscA /„ 

Bnq - - T ^1 - 5-^- - ^1 + -jg-J , (24) 

CN « = -27l 3 -^2rJl 1+ ^2rJ "M^^arJ- (25) 



These coefficients depend on the power of a self-interaction term and its relative size to 
the mass term, which can be calculated numerically [2HI22]- There is a clear distinction 
in the relation between /nl and #nl depending on the contribution from a self-interaction 
term. From the viewpoint of the /nl^nl relation, the model with and without a self- 
interaction term can be regarded as "enhanced <?nl" and "linear g NL " types, respectively. 
Thus once /nl and #nl are determined by observations with some accuracy, it may be 
even possible to probe the form of the potential J5T][22| as well. In Fig. [21 we show the 
case with V(a) = (l/2)m 2 a 2 (pure quadratic potential) and V(a) = {l/2)m 2 T a 2 + \o % jra\ 
(with a self-interaction term). As seen from the figure, we can easily see the difference in 
the /nl~<7nl relation between the cases with and without a self-interaction. 



3.2.2 Pure modulated reheating model 

In th e m odulated reheating scenario [60ll6Tj . the decay rate of the inflaton T fluctuates in 



spac a #10 l because of its dependence on a modulus field a, which acquires quantum fluc- 
tuations during inflation, and large non-Gaussianity can be generated [QUEUES]- In this 
scenario, the level of non-Gaussianity is highly dependent on the assumption of the decay 
rate. Thus we need to specify the dependence of T on a to give a concrete prediction. 
Furthermore, in fact, we also need to specify the form of the inflaton potential and the 
interaction for the decay. But we first give a general expression for the curvature pertur- 
bations ( in this scenario, taking into account the above mentioned uncertainties. £ in 
this scenario can be given by 



r 1 / r r 2 , 



Cmod = A{x)f5cr* + - M(x)-yr + B(x)^j 5a 

+1 (Mx)^ + 3B(x)^^ + C(z)G) Sal (26) 

Here T a = dT/dcr, T aa = d 2 T/dcr 2 and so on. A(x),B(x) and C(x) are functions of 
x = T/H c with H c being the Hubble parameter after several oscillations of the inflaton 
field. Since the number of e-folds after inflation only depends on the quantity x = T/H c , 



* 10 Thus, the reheating temperature after inflation also fluctuates in space in this model, which may 
generate significant isocurvature fluctuations and give severe constraints in some settings [521155] . 
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which can be confirmed by a dimensional analysis of the background evolution equations, 
the coefficients such as A, B and C depend only on x. From the above expression, we 
obtain the non-linearity parameters for this model as 

6 f B(x) i rr aa 

5 /nl - Mxf A(x) VI ' (27) 
54 C(x) 3B(x)TT aa T 2 T aaa 

25 m " Mx) Mx) n + "T3— (28) 

To make some explicit predictions for these parameters, we need to know the values of 
A(x), B(x) and C{x) for a given x, which requires us to specify the interaction responsible 
for the decay and the potential under which the inflaton oscillates. Furthermore, some 
numerical analysis may be needed. However, some approximations are available which can 
give quite accurate results, in particular, for the case of x <C 1 [371165] . 

To evaluate to what extent the non-linearity parameters can be large, here we assume 
that the potential of the inflaton is V oc 2 and the inflaton decays through £ int = —y^tpip. 
Then the coefficients are evaluated as A = —1/6, B = 1/6 and C = —1/3 (for other types 
of potentials and interactions, see [65]). In this case, the non-linearity parameters are 
given by 

j/nl = 6 - 6^, (29) 

|jnl = 3 6 (2-3^ + ^)^ (30) 

To get some numbers for /nl arid <7nl, we further need to assume an explicit form for 
T. Here we take the following form: 

r = r„(i + 4 + /4), (3i) 

where a and (3 are some coefficients, M is some scale and we assume |er| <Mto justify 
the exp ansio n in terms of a/M. With this form of T, the non-linearity parameters are 



given bv r 11 



* n Although here we assume that the curvature perturbation is totally generated from this mechanism, 
in general, fluctuations from the inflaton would also contribute, in particular to the power spectrum. 
In order that the contribution from the inflaton is negligible at linear order, M/M v \ < 0.2ay/e should 
be satisfied with e being the slow-roll parameter for the inflaton, which is obtained by requiring that 
(mod > Cinf with (i n f being the curvature perturbation generated from the inflaton. Another requirement 
comes from the constraint on the tensor-to-scalar ratio r < 0.1, from which we obtain M/M p \ < 0.01a. 
These relations can be satisfied by taking the scale M appropriately. 
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Thus, for example, if we take a = 1/2 and /3 = — 1, then one can have /nl ~ 45 and 
g NL ~ 430. 

In fact, we can find a relation between /nl and (?nl independent of the form of T as far 
as the third derivative T ar7r7 is negligible. After some algebra, one can derive the relation 
between /nl and #nl as [37] 

50 

9nl = 10/nl - y • (33) 

Notice that the signs of /nl and ^nl are the same for large values of /nl, which is different 
from the one in the curvaton model. Thus once the signs of /nl and $nl are determined 
from observations, we can discriminate between the curvaton and the modulated reheating 
models (see also Fig. [2]). In both models, the simple cases predict that #nl is of the same 
order of /nl as given in Eqs. ( 1211) and ( |33|) . which is what we call "linear g^" type. 
Thus, as far as the modulated reheating scenario is concerned, larger g NL may suggest a 
non- negligible third derivative of the decay rate Y aaa . 



3.2.3 Modulated-curvaton model 

In the modulated reheating scenario, density fluctuations originate from those of a modulus 
field. Such a field may also play a role of the curvaton at later time after it "modulates" 
the reheating. In this case, the fluctuations become a mixture of those generated from the 
modulated reheating and the curvaton mechanism although the source of fluctuations is a 
single field a. In fact, this kind of possibility has not been well explored in the literature. 
Therefore, we here discuss this scenario in some detail. 

The curvature perturbation ( in this scenario can be written as the sum of the contri- 
butions from the modulated reheating and the curvaton mechanism: 



c 



Ccur Cmod 
2^dcc 



3(T* 



+ 



81a! 



-9rL + 



i / 2 o 3 \ i i r 

(3rd GC — 4r dcc — 2r dcc J — — 



r 



r 2 



5al 



dec 



°' dec 



1 

36 



2TI 

r 3 



3roT<jcr 
r 2 



+ 



5a: 



(34) 



where Ccur and Cmod are the curvature perturbations generated from the curvaton mecha- 
nism and the modulated reheating, and are given in Eqs. (fT7|) and (126)) . Here we assume 
that the curvaton potential is quadratic and the potential for the inflaton is V oc 2 . The 
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non-linearity parameters are then given by 



5 /NL 



4^dec 

cr* r 



3 



3 



r 



4r 



dec 



at 



(3 - 4r dcc - 2r 



dec J 



(35) 



54 

25^ NL 



-36 



4r, 



dec 



2r 3 , 31^ 



dec 



r 3 



r 2 



9CT 3 



(-18 + r dec + 20rL + 6ry 



(36) 



To discuss the prediction of this model more explicitly, we need to specify the form of the 
decay rate. Here we again adopt the form given by Eq. (!3T|) . Since this model includes 
both cases of the curvaton and the modulated reheating, some limits should be identical 
to those models. However, notice that which part gives a dominant contribution to the 
total curvature perturbation C, can be different order by order. Assuming <r*/M -C 1 and 
rdec *C 1, the condition in which C cur > C mo d at the first order can be given by 



Tdec > Ci — , 

where C\ = a/4 is a constant. In the second order, the corresponding condition is 

r d ec>|C 2 |(^) 2 , 
where = (2/3 — a 2 ) /A. For the third order, it is given by 

CT*\3/2 



-dec > (£) 



(37) 



(39) 



where C3 = (a 3 — 3a/3)/8. As far as the parameters a, {3 are 0(1), the constants C±, C2 and 
C3 are also 0(1). From the inequalities above, one can easily notice that there can exist, 
for example, the case where the curvaton part dominates over the modulated reheating 
part in the second order, however, in the third order, the modulated reheating part gives 
a more contribution. When a, (3 = 0(1), the parameters r dec in the curvaton and a*/M in 
the modulated reheating totally determine the size of £ in each order. From the conditions 
Eqs. ( 137)) . ( 138)) and (139)) . we can divide the parameter space in the <7*/M-rde C plane into 
four regions, which is shown in Fig. [3J Since we consider the case of cr*/M <C 1, for a fixed 
value of cr^/M, the corresponding values of rd ec for each region are: 



^dec,4 < IC2I < 7"dcc,3 < s/\C^\ (^j ' < r d ec,2 < \C\\ < r deCjl 



(40) 



where the subscript % which appears in rd ec ,i indicates the region shown in Fig. |3l Region 1 
corresponds to the case where the curvaton fluctuations always dominate over those from 
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the modulated reheating. In Region 2, only at linear order (^l d > (cm holds, but in the 
second and third orders, Ccur > Cmod an d (cur > Cmod are satisfied. In Region 3, only in 

(2) (2) 

the second order, Qui > C mc , d holds, but in linear and the third orders, fluctuations from 
the modulated reheating give dominant contributions. Region 4 corresponds to the pure 
modulated reheating case. 

10~ 2 
10~ 3 
10" 4 

Uec 10 " 5 

10~ 6 
10~ 7 
10" 8 

10~ 5 10~ 4 10" 3 10" 2 

cr*/M 




Figure 3: Parameter regions where the curvaton fluctuations dominate over those from 
the modulated reheating in all orders (Region 1), and the fluctuations from modulated 
reheating dominate the curvaton ones, in the first and the third orders (Region 2), only 
in the second order (Region 3), and in all orders (Region 4). Note that here "all orders" 
just means "up to the third order" since we consider fluctuations only up to this order. 
We used a = ~ and (3 = — 1 for illustrative purpose. 

Now we look at each case in order. In Region 1, r^cc should satisfy r^cc > \Ci\cr*/M, 
which is required to have ( cur > Cmod at all (up to the third) orders. Here we assume that 
cr* > ^cr* = H*/2ii, which guarantees the linear order perturbation always dominates over 
the second order one. Since Ccur should be equal to the observed amplitude of fluctuation, 
we require £ cur ~ 10~ 5 , which yields another condition r^cc ~ 10~ 5 <7*/i/*. By eliminating 
rdec from the above two conditions, we find that M must satisfy a bound M > 10 5 |Ci|-£/*. 

In Region 2, although the curvature fluctuations from the modulated reheating give 
more contribution than those from the curvaton at linear order, the curvaton dominates 
over the modulated reheating both in the second and third orders. In this case, the relation 
between /nl and #nl can be given by 

SNL ~ 37^/nI 2 - (41) 

Although #nl is proportional to / NL , the size of #nl cannot be very large, comparing to 
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that of /nl, due to the suppression by rd ec - (Notice that r^cc ^ cr*/M in this region.) 
Furthermore, parameter sets in Region 2 generally give too large /nl- The expression for 
/nl in Region 2 is approximately given by 

/nl ^ 5- 77 ■ ( 42 ) 



5 J1 ^ a 2 VM 
r^ec should satisfy the inequality of Eq. ( 140"]) , from which we can find 

6 6 6 

/ . < 7/nL < TTT- (43) 

Since here we assume that cr*/M <C 1, /nl is expected to be large. For example, if we 
take cr*/M < 10~ 3 , /nl ~ 200, which is already outside the current limit. But, in fact, by 
fine-tuning some parameters such as a and f3, we can have the situation where /nl < 100. 
For example, if one takes (3 ~ a 2 /3, this choice of parameters effectively enlarges Region 
2 since C3 ~ (see Eq. ( 140]) ). We show such a case in Fig. [2j 

In Region 3, the relation between /nl and #nl can be written as 

/<t*\ 2 /nl 

£nl ~ 77 • ( 44 ) 



In this region, the inequality (er*/M ) 2 /rd cc < 1 should be satisfied (see Eq. (HQ])), thus 
#nl tends to be smaller than /nl as in the case of Region 2, which can be regarded as 
"suppressed <?nl" type. 

Region 4 corresponds to the limit of the pure modulated reheating case. Thus this 
region follows the argument in the previous subsection 13.2.21 

As a final remark on this scenario, it should be noted that the linear order term can 
be canceled to vanish for some parameter values of r^ cc and a*/M even if a* > 5a* is 
satisfied. In this situation, by introducing the fluctuations from the inflaton, this scenario 
becomes like the ungaussiton model which is discussed later. 



3.2.4 Inhomogeneous end of hybrid inflation 

In this subsection, we discuss a hybrid inflation model, where a light scalar field other than 
inflaton exists and fluctuations of such a field give cosmic density perturbations [BBTTTO"] . In 
usual hybrid inflation models, the inflationary phase ends due to the tachyonic instability 
of a waterfall field and the end of inflation homogeneously occurs. If, however, such a 
waterfall field is coupled with a light scalar field, then the fluctuations of the light scalar 
field can drive the inhomogeneous end of inflation, from which the curvature perturbation 
can be generated. 

Here let us consider one of the simplest models introduced in Ref. [68], where the 
potential for the model is given by 

V = \ (j - X 2 ) 2 + ^ W + + ^/ VV + \ml a 2 , (45) 
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where g, f and A are some coupling constants and v is a vacuum expectation value (VEV). 
77i0 and m a are the masses for the inflaton and a light scalar field a, r espectively. \ 
is the waterfall field. Here we assume that only a acquires fluctuation: 
mass squared of the waterfall field is 



#12 



The effective 



2 i 2i2 i r2 2 

-v +g (f) + f a 



and then a critical value of the inflaton field corresponding to m 2 = is given by 

(Per = • 

9 



(46) 



(47) 



Here, we have assumed that v 2 > f 2 a 2 . Notice that the critical value of the inflaton can 
fluctuate due to the fluctuations of the light field a. Assuming that the inflation ends 
at = cr because of the tachyonic instability of the waterfall field (hence we adopt a 
sudden-end approximation), the total e- folding number during inflation can be estimated 



as 



N 



M 2 y 



V 



(4c 



where we have used the slow-roll approximation. Since the fluctuations of the light field a 
affect the perturbation of the e-folding number, which corresponds to the total curvature 
perturbation (, through the fluctuation of the critical value of the inflaton cr , we have 



dN 



-a 



da 

d 3 N 

Mr 



(T :k 



d 2 N 



da 



ON d 2 



da 2 



da! 



da 



d 2 N 

'Mr 



da 



d 2 



da 2 



dN d 3 <p CI 
d(f> CT da 3 



5at 



(49) 



Hence, the non-linearity parameters are given by 



6 _ 

s /NL " ~m 



1 



b> 2 

cr 



54 
25 



N, 



9nl 



N 3 



N 3 



rcr 

V 2 



1 



rcr 

V 3 



(50) 
(51) 



where the prime denotes the derivative with respect to a. Since in the ordinary hybrid 
inflation model the slow- roll conditions are satisfied until <fi = <f) CT , we find that the first 
terms in Eqs. ( 150]) and (jBTj) are suppressed by the slow-roll parameters. In order to real- 
ize the large non-Gaussianity in this model, large contributions from the second terms in 



# 12 For the case where fluctuations of <j> also contribute to C, which becomes a multi-source case, see the 
discussion in Sec. 13.3.31 
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Eqs. ( 150]) and (!5T|) or the third term in Eq. (l5Tj) are needed. Neglecting the terms sup- 
pressed by slow-roll parameters, the expressions for non-linearity parameters are reduced 
to 

^/nl * -v 7 ^-^ , (52) 

5 0'cr 

54 18 0'" 

T^nl ^ -(2e cr - ??cr)— /nl + 2e cr-f| , (53) 
where e and 77 are the slow-roll parameters defined by 

and we denote the quantities evaluated at = cr by a subscript "cr." 
From Eq. (117)1 . we have 



/v 



(55) 
(56) 



da 3 g A (j)l r \g 2 (pl 

Then, the non-linearity parameters can be written as 



da 2 g 2 <p cr \g 2 (j) 2 

^ s£-(i|. (57) 



54 2 t; 2 

2^ NL=6r ^7v 



^NL = 6?7 2 r -— - , (59) 



where we have used 77 = M P1 v / 2e/</> since y/2e = M^m^/V and 77 = M^m^/V during 
inflation in this model. Furthermore we have also adopted the approximation e cr <C r] cr <C 
1, which is satisfied as in a usual hybrid inflation. Hence, when v 2 ^> f 2 cr 2 /ri cr , large /nl 
can be generated in this scenario. Notice that, in general, fluctuations from the inflaton 
would also exist, which gives another condition for neglecting the inflaton fluctuations: 
f 2 a > g 2 (J) cr (ecr/e*) 1 / 2 . This relation can hold by taking the value of a appropriately. 
Then, we can find the following relation between / NL and <? NL : 

#nl = ^ry /nl • (60) 

Compared to the modulated reheating scenario, the ratio of the non-linearity parameters 
S'nl/Znl is suppressed by the slow- roll parameter 77 in this scenario and thus this model 
can be regarded as "suppressed (7nl" type. Another simple model introduced by Alabidi 
and Lyth [70] also predicts that <7nl//nl is suppressed by the slow- roll parameter. 
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3.2.5 Inhomogeneous end of thermal inflation 



Primordial fluctuations can also be generated by modulating the end of thermal inflation 
[7T] (see also [72] for the discussion on the generation of primordial fluctuation from 
inhomogeneous cosmological phase transition). Here we briefly review the mechanism 
following [71] . 

The thermal inflation [73] can be realized by using a flaton field <p, a flat direction in 
supersymmetric theories, whose effective potential can be given by 



(gT 2 



m 2 ) 



A 1 



6M^ 



(61) 



where (l/2)gT 2 is a thermal correction to the potential due to the interaction of the flaton 
field with some particles in the thermal bath with g being an effective coupling between 
them and T being the cosmic temperature. When ^fg T > m, the flaton field is trapped 
at the minimum = 0, then at some time when the energy density of the flaton given 
by Kons becomes larger than the background energy density, the Universe is dominated 
by the false vacuum in the potential, which drives a mini-inflation. However, when the 
temperature decreases down to T c = m/y/g, the flaton rolls down to the VEV, then the 
mini-inflation ends. The total number of e-folding from the time ij when the mini-inflation 
begins to the time tf when the flaton rolls down to the VEV is given by 



AT(t c ,t i: 



Hdt = In 



tin 



-In 



T 
T~ 

-L ir 



(62) 



where T[ n is the temperature when the mini-inflation begins. If the coupling g depends on 
a light scalar field a as g = g{<r), the coupling g can fluctuate since a light scalar field a 
can acquire quantum fluctuations during inflation, which gives rise to fluctuation in the 
number of e-folding, that is, the curvature fluctuation. By utilizing the 5N formalism, the 
curvature perturbation can be given by 



( = 5N: 
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from which we can derive the non-linearity parameters in this model as 



6 
5 
54 
25 



/] 



NL 



3NL 



- 1 



g"/g 

1(9'/ 9? 

g" /g _ 3 g"/g 
(g'/gf (g'/gf 



5*1 , (63) 



(64) 
(65) 
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where the prime represents the derivative with respect to cr*. When g is negligible, the 
following relation holds between / NL and #nl : 

9nl = -y/NL - — • (66) 

Thus this model can be considered as "linear g NL " type in terms of the /nl~<7nl relation. 
Although the mechanism is quite similar to the modulated reheating scenario and the 
inhomogeneous end of hybrid inflation discussed in the previous subsection, the relation 
between /nl and #nl is different. In particular, the relative sign between /nl and c/nl 
differs from the one in the modulated reheating although the size of /nl and #nl is almost 
the same order. On the other hand, the size of #nl is larger than that for the case of 
the inhomogeneous end of hybrid inflation, where <?nl//nl is suppressed by the slow-roll 
parameter. 



3.2.6 Modulated trapping mechanism 

During inflation, the resonant particle production can occur because of the coupling of the 
inflaton to other particles [7H476] . If the relevant coupling A depends on a light scalar field 
a, dubbed as modulaton in [77], the curvature perturbation can be generated through 
fluctuations of a which originate from quantum fluctuations during inflation. Such a 
mechanism is called modulated trapping, which was proposed in [77j . 

Here we briefly discuss this mechanism, following closely [77], and give some relation 
between /nl and $nl in the model. Let us consider the following interaction for the inflaton 
0: 

Ant = ~N"(m - \<P)xx, (67) 

where x is a fermion coupled to the inflaton and M is the number of species of particles 
with the same mass. When the inflaton crosses the value <fi = <fr pp = m/X, x held becomes 
effectively massless, then they are resonantly produced. At this moment, the occupation 
number for x abruptly increases from zero to 

\3/2 

^ PP = ^|0pp| 3/2 , (68) 

where "pp" indicates that the quantities are evaluated at the time of particle produc- 
tion. Once the particles are produced, its number density just decreases with the cosmic 
expansion, thus after its production, the number density can be written as 

n(t) = n vp (— y ) e(t-t P p). (69) 
V a pp/ 

By using the above expression and adopting the Hartree approximation in the equation 
of motion for 0, we can evaluate the effect of the backreaction of the particle production 
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to <fi. The equation of motion for can be written as 



+ m4>+ d ^jp-=M\n p J— ) 3 Q(t-t pp ). (70) 



Assuming that the particle production and its dilution occur in a short time scale compared 
to the cosmic expansion during inflation, and dV (cj)*) / d(j) can be regarded as constants 
during the particle production. Defining 

A0(t)=0(t,A^O)-0(t,A = O), (71) 

then, from the equation of motion above, we obtain 

A0(t > t pp ) = Af\n pp ex P [-3H pp (t - t pp )](t - t pp ). (72) 

Now we evaluate the curvature perturbation by using the SN formalism. The number 
of e-folding which is attributed to the particle production can be given as 

AA r(A^o) = _# pp ^. (73) 
0pp 

By integrating A0 given in Eq. (1721) from the time of particle production, we obtain A</> 
as 

A0 = jf A<pdt = ^* (74) 
Thus AN^ ^ can be evaluated as 

Ai ^»> = -U.^L = ^MH, (75) 
0pp 18tt s H pp 

from which the curvature perturbation due to the fluctuations of a is given by the formula 
in the 5N formalism: 

C = (AN^)ja* + 1 (AiV^o)) ^ ( ^ )2 + 1 (AiV(^)) _ (5a,) 3 . (76) 

To give some explicit expressions for the non-linearity parameters such as /nl and (/nl, 
we need to assume the functional form or a dependence of A and m, which is somewhat 
model-dependent. As discussed in [77], in the case where A is independent of a, but the 
mass depends on a as m = ga, the non-linearity parameters cannot be much larger than 
unity. Thus we do not consider such a case here. When the coupling A and the mass m 
depend on a as A = a/M and m = ga with M being some energy scale, /nl and #nl are 
evaluated as 

6 9 

t/nl = — (77) 
5 bep 

54 27 

25^ NL " 25^' (78) 
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where f3 is the "efficiency factor" , which is defined and given as 



= 




(79) 



Since the particle production occurs at the cost of reducing the kinetic energy of the 
inflaton, cannot exceed 1. If we take (3 ~ 0.01, /nl can be as large as /nl — 55 [77]. 
Furthermore, from the above expression, we can find the relation between /nl and ^nl as 



which is "enhanced <7nl" type. In fact, the relation g NL ~ f^ L holds even if we assume 
somewhat general type of functional form for A = (a/M) p . Thus this mechanism predicts a 
large values for #nl relative to /nl compared to some other "modulated coupling" scenarios 
such as modulated reheating and inhomogeneous end of thermal inflation, which give 



3.3 Multi-source model 

Although the origin of density fluctuations is usually assumed to be a single source, in 
general, it is possible that multiple sources can be simultaneously responsible for density 
fluctuations [5T | |65 |[78H8T] . For example, in the curvaton scenario, the curvaton field alone 
is usually supposed to generate the curvature perturbation, however, even in this scenario, 
the inflaton should exist to drive the superluminal expansion at the early epoch and it can 
acquire quantum fluctuations. Thus in general, both the inflaton and the curvaton fields 
can be responsible for cosmic density fluctuations today. In a case where fluctuations 
from multiple sources can simultaneously give a sizable contribution to observed ones, 
the relations among the non-linearity parameters are different from those discussed in 
Section 13.21 

Here we discuss the case where two scalar fields and a generate density fluctuations. 
In this case, the curvature perturbation can be generally written as 



3NL = 





9nl\ ~ |/nl 



C = N+Sfa + N a 8a* 
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Then the non-linearity parameters can be given by 
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R represents the ratio of the power spectrum for the curvature perturbation from fluctu- 
ations of 6 and a: 



Pf\k ief ) 
P^ ( k rc f ) 



R 



(86) 



where 



P^\k vci ) = N*P s (k Tei ), P^(k vci ) = N 2 a P s (k Ici ), (87) 

with Ps being defined in Eq. (jSJ) and the power spectra are evaluated at some reference 
scale fcref- Then, we obtain the following relation, 



tnl 



where the ratio R is defined as 



1 + iA /6 , 

IT) U /nl 
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These formulae can be applied for a general two-field case. 
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3.3.1 Mixed model with inflaton fluctuations 



In a special case where </> is the inflaton field and a is some other light field, (p does 
not contribute to non-Gaussianity because non-linearity parameters are suppressed by the 
slow-roll of <f>. Then, the curvature perturbation can be written as 

(4>) , A°) 



c = c w + c 



I o 



(91) 



where and C^ a \ represent the contributions to the curvature perturbation from Scj) and 



5a, respectivel y^ 13 ! . In most cases, the inflaton and the other light scalar field can be con- 



sidered to be uncorrelated and non-linearity in the inflaton sector is generally suppressed 
by the slow-roll, here we set fjf^ = 
Furthermore, in this case, we can a 
relation among /nl, tnl and #nl as 



by the slow-roll, here we set fjf^ = f^f^ = 0, which indicates that -R/ NL = and R = R. 
Furthermore, in this case, we can also set = g^L^ = 9nl^ = 0- Then we have the 
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TNL=(^) (I/nl) , (92) 



and 
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(93) 



This equation can be applied for a scenario where one of two sources of fluctuations alone 
contributes to non-Gaussianity. Also notice that when R is small, it corresponds to the 
situation where tnl 3> (6/nl/5) 2 . Thus large tnl/ /nl ma y indicate that (almost) Gaussian 
field gives a lar ge c ontribution to the power spectrum, but non-Gaussianity comes from 



the other source* 14 



In most models, the relation between and g^l for a single field sector can be given 
by the following form: 

g^^C.JfSY + C^, (94) 



where C*i CT and Cio are numerical constants and p is the power of f^. For example, 
in the curvaton model with a quadratic potential, these parameters are given as C\ a = 
— (10/3), C 2a = —(575/108) and p = 1 as can be read off from Eq. (j2T|) . By assuming the 



* 13 In some literature, the notation ^ indicates that the curvature perturbation on the slice where the 
energy density of (f> is uniform. Here £^ just represents the contribution to £ from Scj). 

* 14 Here it should also be noted that, for certain parameter values of this kind of mixed models, the 
second order fluctuations give the leading contribution in In this case, the model becomes like 

"ungaussiton" model, categorized as "constrained multi-field type," which will be discussed in the later 
section. 
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form of Eq. (194)) for the a sector, we can write the relation between /nl and <?nl for the 
mixed model with the inflaton as, 

R \ 3 ' 2 " .„ _ ( R 



^ = c >°{ttr) ^ +c ^{ttr) ■ (95) 

We can also write this equation by using tnl and eliminating R as, 



fi \ 6-4p p6-3p /«\ 6 f6 

c ^ + (j>) c *k- (96) 



As an explicit example, first we give the formula for the case of the mixed curvaton and 
inflaton fluctuations. As mentioned above, even in the curvaton scenario, the fluctuations 
from the inflaton can also give some contribution to the curvature perturbations as well 
as those from the curvaton. Here we denote a as the curvaton field. As discussed in 
Section 13.2. H assuming that the potential of the curvaton is quadratic, and g^l are 
related as Eq. (121]) for r^cc <C 1. Putting the relation into Eq. ( 193]) or ( 196]) . we find the 
relation among the non-linearity parameters for the case with r^cc < 1 as [78] . 

24 /£ L 9936 /£ L 

s'nl = — z — — 3-. (y<j 

5 r NL 625 rf^ L 

This relation can be rewritten as 



10 / R \ 575 / R 

^ NL — 3~ v r+~R ) — 108 v r+~R ' {9i 



which recovers the pure curvaton case Eq. (12 lj) in the limit of R — > 00. When R ~ 1, $nl 
is of the same order of /nl, but when i? -C 1, g NL is suppressed by the factor R. 

A similar relation can be obtained for a mixed model of the modulated reheating 
and the inflaton fluctuations. As discussed in Section 13.2.21 when the term with Y aaa is 
negligible, is related to g^l as in Eq. ( 133]) . By adopting the relation, we obtain the 
corresponding equation in this model as [65] 

"^.i^^of * * V, (99) 

5 r NL 625 r^ L V 1 + ^/ 3 \1 + R J 



which reduces to the pure modulated reheating case (1331) in the limit of R — > 00. Notice 
that g NL is again suppressed by R when R -C 1. 

Another example is the modulated trapping model with a sizable contribution from 
the inflaton fluctuations. Here a is assumed to be modulaton. For some functional forms 
of A and m, the relation between and g^l can be given as Eq. (150"]) . In this model, 
the following relation holds: 

2 fl + R\ f2 25 , 1nnN 
9 l-H/m.= 162^,. (100) 
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In general, the inflaton fluctuations can also contribute to the primordial fluctuations 
even if we consider another source. The discussion above is applicable for other models 
regarding a as the one in other single-source models. If the corresponding single-source 
model predicts the /nl-^nl relation of the linear type, the counterpart of its mixed model 
with the inflaton gives "suppressed #nl" type due to the ratio R when i? C 1, although 
|/nl| ~ |#nl| when R~ 1. 

As a final remark, we comment on the spectral index n s and tensor-to-scalar ratio r in 
this type of mixed models. n s and r in this scenario are generally given by 

n, - 1 = -26 - ^f, (101) 
166 (102) 



where e and 77 are slow-roll parameters for the inflaton defined as in Eq. ( 154|) . The limit of 
R — > corresponds to the case where only inflaton fluctuations are responsible for today's 
density fluctuations. 

On the other hand, the limit of R — > 00 corresponds to a single-source model. As 
seen from the above equation, the tensor-to-scalar ratio r becomes very small in this limit. 
Thus r is considered to be generally very small in models with large non-Gaussianity. 
However, when R ~ 1, large r is possible while non-Gaussianity can also be large. Thus a 
mixed model discussed in this section is interesting in this respect as well. 



3.3.2 Multi-curvaton model 

In a usual curvaton model, there exists only one curvaton field. However, it may also be 
possible that multiple fields can play a role of the curvaton. Some authors have already 
considered this kind of the scenario [82l[83], which is called "multi-curvaton" model. Here 
we discuss this model focusing on its non-linearity parameters. 

In the following, we basically follow the arguments in [83]. Here we adopt the sudden 
decay approximation and assume that two curvaton fields, denoted as a and b, are respon- 
sible for the present-day density fluctuations. We also assume that the curvaton field a 
decays first (i.e., the decay rates for a and b are assumed to be T a > T b ). At an infinites- 
imal time before and after the first curvaton decay, the total curvature perturbation d 
at the first curvaton decay, the curvature perturbations on the constant curvaton a and b 
hypersurfaces ( a and ( b are related as 

(1 - O al - Q bl )e^ + Q al e 3( ^ + Q bl e 3 ^-^ = 1, (103) 
(1 - n bi y^-^ + ftwe 3 ^ = 1, (104) 

where Q a i = Pai / Pi and Q b i = Pbi / Pi with pi , p a \ and p b \ being energy densities of the 
total component, the curvaton a and b at the first curvaton decay, respectively. At the 
time just before the second curvaton decay, the following equation holds: 

(1 - tt b2 )e 4 ^-^ + n b2 e^-^ = 1, (105) 
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where Qb2 = Pvij ' P% with p 2 and pb2 being energy densities of the total component and the 
curvaton b at the second curvaton decay, respectively. Here £ 71 is the curvature perturba- 
tion on the constant radiation hypersurface at the first curvaton decay. (2 is the curvature 
perturbation after the second curvaton decay. Here the subscripts 1 and 2 indicate that 
the quantities are the ones at the time of the first and the second curvatons (the curvatons 
a and b) decays, respectively. 

Since we are interested in the final curvature perturbation after the second curvaton 
decays, we evaluate (2 up to third order, then find non-linearity parameters in the model. 
The curvature perturbations ( a and (b are related to the field perturbations of the curvatons 
a and b, denoted as 5a and 5b, by 

1, . . 2 5a 1 f 5aV 2 f 5a\ 3 . . 

C„ = -io g (i + <y = — +- — , (106) 



3 3a 3 \ a J 9 \ a 

1 2 5b 1 {5b\ 2 2 ■ 



a = _ log(1+4) = ____^_j + , m 

where 5i = 5pi/pi and we truncated at the third order. Here we have assumed that the po- 
tentials for the curvatons are quadratic. Hence the field values a, b and their perturbations 
5a, 5b are regarded as the ones evaluated at the time of horizon crossing. 

By using an iterative method, we can express (2 as the series in ( a and (b ■ Denoting 
the first order parts for the curvatons a and b in Eqs. (11061) and (11071) as £ a (i) and (b(i), 
respectively, we can generally write (2 as 

C2 = C a Ca(l) + C&C&(1) + CaaCa(l) + C&&Cb(l) + Ca&Ca(l)C&(l) 

+Caaa(a{l) + CbbbCb(l) + C aa b(a(l)(b(l) + C ab bCa(l) C&(1) • (108) 

Here the coefficients such as C a and so on are functions of f a i, fbi and fb2 which are defined 

as 

qo . 

(109) 
(110) 
(111) 



fal — 


3£l a i 


4 — ^ol — ^61 


fbl = 


3£lbi 


4 — fi a l — ^61 


fb2 = 


3^62 


4-fi fe2 ' 



These quantities roughly represent the ratio of the energy density of the curvatons a and 
b to the total one at the first and second curvaton decays. 

Since the expression for the coefficients are very complicated in general, here we con- 
sider only some limiting cases where non-Gaussianity can be large. We give the full 
expression of the curvature perturbation in Appendix [A] 
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• Case I: Both curvatons are subdominant at their decays 

When the curvatons are both subdominant at their decays, which corresponds to the 
case of fai, fb2 <?C 1 (and fbi < jm being implicitly assumed), the curvature perturbation 
after the second curvaton decay can be given by 

C2 = /alCa(l) + /&2C&(1) + -^Co(l) + ~T~(b(l) 

^"Ca(l) ~~ -^falfb2(a(l)(b(l) ~ ^falfba(a{l)(b(l)i ( 112 ) 

where we have kept f a i and /&2 only at the leading order. 
Then we can find the non-linearity parameters for this 

6 f . , 3 fl + tf 2 K* 

5 /NL -2(^ 1 + /fe W {1U} 



^9 /oi + fb2-K 3 

H (114) 



where is the ratio of the power spectra of the curvatons a and b, denoted as P^ a and 
P^ b , respectively: 

K=& (116) 
With this definition, the ratio R defined by Eq. (1861) in Section 13.31 can be written as 



fai 

In the language of Section 13.31 we have 



R ~ ^-K. (117) 



Regarding the relation between #nl and /nl, we find 



/ Q 3 i + j|g 



:n8) 



y/NL < ~9nl < 10/nl- (119) 

Thus, #nl is of the same order of /nl with the opposite sign in this case and we could 
write as 

2NL = C mc (/ Ql ,/ b2 ,K)/ NL , (120) 
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Figure 4: Contours of — C mc in the f a i~fb2 plane. Here we take K = 1. Notice that, in 
the lower right and upper left region, C mc — > —(10/3), which corresponds to the single 
curvaton case. 
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where C mc is a negative coefficient of 0(1), which slightly depends on the parameters 
fail fa an d K. In Fig. HJ contours of — C mc are plotted in the fa- fa plane, from which 
we can see that C mc is generally 0(1). 

We also find the following inequalities for three non-linearity parameters: 

_ 125 < /nlSnl ^ (m) 
27 r NL 

£ _| (122) 

which may be useful for discriminating this case from other ones. 



• Case II: Both curvatons are dominant at their decay 

It was first noted in [S3] that large non-Gaussianity can be generated even if the 
curvatons are dominant when it decays. This is possible only when there are two curvatons. 
For this case, we assume that fa = 1 and fa <C 1. Although we also assume that fa ~ 1, 
we keep fa up to the leading order, i.e., we expand the expressions for ( 2 and non-linearity 
parameters around fa = 1. The limit K — > oo corresponds to the case where the first 
curvaton effectively homogeneous and only the second curvaton fluctuates. This situation 
is exactly the same as the standard (one-field) curvaton case. However, the case of K — > 
is quite different from the standard curvaton model, which we are going to consider here. 
Taking the limit K — > 0, which corresponds to the case where the second curvaton is 
effectively homogeneous, the curvature perturbation is given by 

5 

C2 = (1 — fa)(a(l) + fa(b(l) + -T7Z TTP - fa)(a(l) + fa(b(l)} 

4(1 - Jb2> 

5 

+ 7o7l TVI^ 1 ~ fb2)Ca(l) + fa(b(i)} 3 , (123) 

t^U - Jb2) 

from which we can find that this case effectively reduces to the single-source case by 
regarding (1 — fa)C a (i) + faCb(i) being a Gaussian field. The non-linearity parameters for 
this case can be evaluated as 

6 /NL = 0/1 5 ; r , (124) 



5 J ^ 2(1-/, 



1/2, 
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25 5NL = 2(1-/, (125) 
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Then we have the following consistency relations: 



36 

TNL = —fL, (126) 



25' 

9nl = ^fL- (127) 
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In this case, #nl is proportional to f^ L and hence may become relatively large, which could 
be categorized as "enhanced <?nl" type. 

3.3.3 Multi-field inflation model 

As is well known, the non-linearity parameters in the standard single field inflation models 
are suppressed by the slow-roll parameters evaluated at the time of horizon crossing. Due 
to the special property of the single field model that the curvature perturbation is conserved 
on super-horizon scales, the breakdown of the slow-roll conditions at the end of inflation 
does not induce additional corrections to the non-linearity parameters. One may then 
expect that the situation changes either if we go on to multi-field inflation models where 
more than one field contribute to the inflation dynamics, or if more than o ne fi eld enter 
the game at the end of inflation where the slow- roll conditions are violatec P^ . This is 
what we want to address in this subsection. 

We first briefly review the general argument which shows that, from a naive order 
counting, /nl, Tnl and <?nl generated during the slow-roll inflation are suppressed by the 
slow-roll parameters. Then we consider a counter example given in [83H86] , where the non- 
linearity parameters can be large during the slow-roll inflation. Then, we also consider 
another multi-field inflation model called multi-brid inflation model, where large non- 
Gaussianity is generated at the end of inflation. 

We set Mpi to be unity only in this subsection to simplify the equations. 



General slow-roll multi-field inflation model 



Here we just present the formula for /nl in a (multi-field) slow-roll inflation model 
with canonical kinetic terms. The formula for /nl at some time tf during inflation when 
the slow-roll approximation is valid can be written in terms of the potential of scalar fields 
V as ETMEl 



5 



iV^6 a (iV f )e b (iV f ) + / dN N a (N)Q a bc (N)Q b (N)Q c (N) 



A'., 



:i28i 



where N" = N a (N*) and N^ b is given by 



f fv 2 



N=N { 



Vab 



V c V d V a V b , V a V b A Vc(aV b) V< 



(v 2 ) 



' 2 \ 2 



-Vcd + 



V 2 



V 



(129) 



N=N { 



Here we introduced a notation for tensor indices as = {l/2){t ao +t ba ) and V a = dV/d(f> a , 
Vab = Vab/V and V' 2 = V a V a . All the quantities on the right hand side in Eq. (1 1 2 8 [) are 
evaluated on the background trajectory in field space and the e-folding number is used as 



# 15 The special cases belonging to this class have already been discussed in 13.2.41 and !3.2.5l 
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the time coordinate (iVf and TV* are respectively the e-folding numbers corresponding to t{ 
and i*). Here, we have taken the final hypersurface at N = Nf to be the V = constant one, 
which is approximately corresponding to the uniform energy density one under the slow- 
roll approximation. The index is raised by the inverse of the field space metric, which is 
assumed to be S ab here. To define N a and a , it is convenient to introduce the propagator 



AW N') 



N 



Texp / P(N")dN' 



N' 



(130) 



where T means that the matrices P a b are ordered in time when the exponential is expanded 
in power of P a b and 

pvjv) = -y + -^. (isi) 

Using this propagator, N a and 6 a are defined as 

N a (N) = N { b A b a (N { , N), Q a (N) = A\(N, N>)N b „ (132) 

with 



N, 



(133) 



b V V a 

v a V N=N{ 

N a (N) and a (A^) respectively satisfy the following equation of motion: 

^N a (N) = -P a b N b (N), (134) 

^e a (N) = p b a e b (N). (135) 

The boundary conditions are given by N a (Nf) = N l a and 9 a (iV,) = N a (N*), respectively. 
The three point interaction Q a bc is given by 

V a h V a uV r V a Vh V a V hr V a V h Vr 

Qbe(N) = ^ + -^ + -^ + -^-2—^. (136) 

In a naive sense, the slow-roll conditions require that the potential of the inflaton is a 
smooth fu nctio n of <p a and hence higher order differentiations with respect to <p a are more 
suppresse d^ 16 . Using the similar slow-roll parameters as defined in Eq. ( )54|) and assuming 
those for every field as of the same order 0(e), P\ and Q a bc can be estimated to be of 
0(e) and 0(e 3 / 2 ), respectively. Since the duration of the inflation is roughly estimated as 

w - (afe)-°(g)- ^)- ^ 



* 16 The assumption of this kind of scaling would be natural to obtain an almost scale-invariant spectrum 
although it is not strictly required [55] , 
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the exponential factor in A a b defined in Eq. ( 11 30[) is 0(1). Furthermore, the value of 
is estimated as 



v b v h 



0(e 



(138) 



Hence, one can see that N a and a are 0(e" 1 ^ 2 ). Substituting these estimates into 
Eq. (I128p . we find that /nl = 0(e). This rough estimate indicates that /nl is typically 
smaller than 0(1) within the range of validity of our present approximation. 

However, when there are some large hierarchy between the slow-roll parameters, which 
is possible for a multi-field case, /nl can be much larger than 0(e). For example, let 
us consider the case where there is huge difference in size between e\ and e 2 (which are 
slow-roll parameters for <pi and <f> 2 , respectively). Assuming ex 3> e 2 , /nl can be estimated 
as /nl ~ ^(e) x (^1/^2), which can be large even though e l5 e 2 < 0(1). Hence, a multi-field 
inflation model can generate large non-Gaussianity although a typical estimate of /nl is 
0{e). 

Furthermore, the formulation given above is applicable only when the slow-roll condi- 
tion is satisfied. However, when the slow-roll condition is violated, in particular, at the 
end of hybrid inflation, one could have large /nl- Some explicit example models of large 
/nl will be discussed in the following. 

The above discussion on /nl, of course, can be applied to the higher order non-linearity 
parameters, i.e., tnl and #nl- Following Ref. |3lJ, t~nl and #nl are respectively given by 

Tnl = M a (N*)n a (N«) , (139) 
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fi'NL 



1 



N abc®( N f)® b ( N f)® C ( N f) 
Nf 



dNN a (N)S a bcd (N)G b (N)G c (N)G d (N) 



As 



+3 / dNQ a (N)Q a bc (N)e\N)e c (N) 



where Q a (N) is a new vector variable obtained by solving 

^a(N) = -n b (N)P b a (N) - N b (N)Q b ac (N)Q c (N) 
with the boundary condition Q a (Nf) = N* b Q b (Nf). N { abc is given by 



(140) 



(141) 
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a VaVbVc VaV b V c V d V e V f V d ef , v d v e v a v b v c 



V 3 



(v 2 ) 



,2\3 



V "'^V (V' 2 ) 2 

de VdV e (aV b V c ) 



(v 2 ) 3 

VbVc) 



V 



(v 2 ) 



(v 2 ) 



,2\ 2 



y d v {a v bc)d v b v d c) v d v e v {a 



V 



,2 



V 



- QVd(a- 



V 



, 2 +QVde- T . T//2 



w 



V c 



f 

be) 



(142) 
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and the four point interaction S a bcd (N) is given by differentiating Q a bc with respect to 
cf) a . Then, based on the standard slow-roll approximation, where V a b/V = 0(e), V a b c /V = 
0(e 3/2 ) and V abcd /V = 0(e 2 ), the orders of the non-linearity parameters tnl and #nl are 
respectively estimated as 

tnl = 0(e 2 ) , <7nl = 0(e 2 ) . (143) 

From this rough estimation, we find that tnl and #nl are more suppressed by the slow-roll 
parameters than / NL . 

As in the case of /nl, the magnitudes of tnl and #nl are not necessarily suppressed like 
0(e 2 ) and can be much larger than the rough estimation in some situations. Furthermore, if 
slow-roll conditions are violated, we may have a chance to generate large non-Gaussianity, 
which is not accommodated in this formalism. Indeed, some explicit multi-field inflation 
models have been constructed that can produce large non-Gaussianity. In the following, 
we consider the generation of non-Gaussianity in those models. 



Two-field slow-roll inflation 



In Ref. [84H86] , the authors have shown that the large non-Gaussianity can be generated 
even during the slow-roll inflation. As as example, we consider the potential of the form: 

V((f>!, <t> 2 ) = V M exp Qr/!0 2 + ^ 2 2 ^ . (144) 

In the following, we assume that |?7i0 2 | <C 1, 1 772^2 1 ^ 1 and VJ n f is constant. In fact, the 
above type of potential can lead to large negative /nl in order to have a red-tilted spectral 
index, which is required to be consistent with current observations. Since positive /nl is 
favored at this moment, this model might not be so attractive in this respect, however, we 
briefly discuss this model as an example, in which large non-Gaussianity may be generated 
even during a slow- roll inflation when two fields exist. 

By using slow-roll solutions, we can write the inflaton field values during inflation as 

0i = <pue~^ N , <\> 2 = <\> 2 ^ N . (145) 

Here, we consider only the slowly-rolling phase and do not specify how the end of inflation 
is triggered. Hence we evaluate the curvature perturbation on the uniform energy density 
hypersurface, which can be determined by V = constant hypersurface as in the "general 
slow-roll multi-field inflation" model discussed above. 

The spectral index and the tensor-to-scalar ratio are calculated as 

fa _ 2ee 2, lj V )6ie - 2 , 1 7V + ( _ 2ee 2 V2 N )€2e -2 V2 N 

Hs ~ 1 = ~ 2e * + 2 2e ie ~^N + 2e2e -2 V2 N . (146) 

(147) 



2e 1 e- 2 n N + 2e 2 e~ 2 ^ N ' 
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The non-linearity parameter /nl can be evaluated as 



6 , -e (rjxex + r} 2 e 2 e^-^ N ) + (2/e)e 1 e 2 (r ]l e 2 + ^(l - e 2 ^~^ N ) 2 

5 /NL " (e 1 + e 2 e^-^N ? > ( 148 ) 

where the slow-roll parameters e\ and e2 are defined for <pi and 02, respectively, in the 
same manner as in Eq. (|54p and e = ex + e 2 . It has be en s hown that large non-Gaussianity 



can be generated when the following condition is met r 17 



2 

< 1. (149) 



4? + 0| ei + e 2 



In this case, we can have simpler formulae for n s , r and /nl as follows 



2(771 + Rm) 

n s - 1 ^ W ^ + jR / 2J , (150) 

166* 

151 



1 + iT 



5 /NL " (1 + i?) 2 ^ 26 ' (152) 

where i? is defined as R = N^JN^ = (e 2 /ei) e^-^^see Eq. (p])). From the above 
expressions, we find that 7/2 > should be satisfied to have large positive /nl- However, 
large /nl requires rji > r/ 2 to have a big factor from e 2 ( m ~ m ^ N , which gives bl ue-t ilted 
spectral index. Thus, to have a consistent value of n s with current observation s^ 18 ! , /nl 



should be negative even though its size can be very large. Thus, in this sense, this model 
might not be a promising one from the viewpoint of current observations. However, one 
could consider another potential, which may give different predictions for n s ,r and /nl- 

As a final remark, we mention the relation between /nl and #nl, which holds for the 
case considered here |85j : 

= y yTr /nl - (153) 

Thus, this model can be regarded as "suppressed </nl" type. For tnl, we do not find 
anything more than the relation given by Eq. (19211 . 



* 17 In fact, there is another case where large non-Gaussianity can be generated. However, such another 
case follows the same argument here by changing some parameters appropriately because of the symmetry 
of the potential [85] . 

# 18 Current limit for the spectral index is n s = 0.963 ± 0.012 (68%CL) pp. 
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• Multi-field hybrid (Multi-brid) inflation 

Here, we discuss a multi-field hybrid inflation (dubbed as multi-brid inflation [90] ) 
model. The analysis of the non-Gaussianity in this model is similar to that in the previous 
multi-field slow-roll inflation model, where multiple (inflaton) fields can affect the dynamics 
of inflation. The difference between this multi-brid model and the previous multi- field slow- 
roll model is mainly the assumption on how the inflation ends: in the multi-brid model, 



the end of inflation is characterized by an ellipse in the field spac a^ 19 ! . In the following, 
we discuss two types of potential in the multi-brid inflation model. 



A. Quadratic potential model 

Let us first consider the model whose potential is approximately quadratic [93] : 

V = V exp (Irntl + \r) % <& ] , (154) 



2 ,T1 2 



where 



K) = ^G(0 1 ,0 2 ) X 2 + ^(x 2 -y) , (155) 

^(01, 02) = <?i(0i cos a + 02 sin a) 2 + g 2 {— 0i sin a + 02 cos a) 2 . (156) 

Here 7/1 and 772 can be regarded as the masses squared for the inflatons 0i and 02 
normalized by Vo, respectively. \ is a water-fall field and its VEV is characterized 
by v. The angle a corresponds to the rotation of the ellipse of G(0i,0 2 ) = const, 
in the field space relative to the 0i axis. The slow-roll equation of motion for 0j is 
given by 

^ = rn<t>i, (157) 
where N is the number of e-folds and dN = —Hdt, from which we obtain 

N = — In 0i In 0i f . (158) 

The inflation ends when the following condition is satisfied, 



v 



- ^1 (0i,/ cos a + 2) / sin a) 2 + g 2 (—4>ij sin a + 2 j cos a) 2 , (159) 



* 19 There are several related studies [9TJ[92]. This model, in fact, includes the inhomogeneous end of 
hybrid inflation model, discussed in section 13.2.41 as a single source limit. 



36 



where faj denotes the field value at the end of inflation. In order to parameterize 
(pi j, we introduce an angular parameter 7 defined as 



9i 

v 

92 



cos 7 = <pij cos a + <p 2 j sm a, 
sm 7 = —0i ,/ sin a + 2> / cos a. 



These equations imply 



v v 

n r = — cos a cos 7 sm a sm 7, 

9i 92 

, V . V 

!>2 f = — sm a cos 7 H cos a sm 7. 

Pi 92 



(160) 
(161) 

(162) 
(163) 



The angle 7 represents the position of the inflationary trajectory at the end of 
inflatio n^ 20 ! . From the equation of motion, we can show that is a constant 

of motion, from which we can derive the following equation: 

— In 0i In 02 = — In 0i f In 02 f 

Vi V2 Vi V2 



iln 

m 



9x92 



(<72 cos a cos 7 — #i sin a sin 7) 



-iln 



9\92 



(g 2 sin a cos 7 + <?i cos a sin 7) 



. (164) 



By expanding N given in Eq. ( 11581) . we can obtain 57V, up to the third order, as 



5N = — 



Vi L 0i 



501 



Id 2 In 0i, 



<9 7 2 



(5(1)7 + 5(2)7)' 



d In 0i 

1 <9 3 In 1) 
6 c?7 3 



- (<*(i)7 + ^(2)7 + £(3)7) 



(5(i)7)' 



(165) 



where 5(i)7, 5(2)7 an d 5(3)7 are respectively the first, the second and the third order 
perturbations of 7. Since they can be described in terms of 50; by the relation given 
by Eq. f)164p . the resultant expression of 5N or the curvature perturbation ( with 



#20 rp^g definitions of a and 7 are shown schematically in Fig. 1 of 
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respect to 5< 
C = SN 



is given by 
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where a prime denotes the derivative with respect to 7 , and Ai,A 2 ,F and G are 
defined as 
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77 2 Ai - 771A2, G(j) = A[A 2 - AxA' 2 . 
The spectral index and the tensor-to-scalar ratio are respectively given by 
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Here it is worth noting that the tensor-to-scalar ratio r can be relatively large even 
when /nl is large in this model [93] • As mentioned in Section 13.3.1} the tensor-to- 
scalar ratio r tends to be very small in most models generating large non-Gaussianity. 
However, the multi-brid inflation model is one of those which can realize large tensor- 
to-scalar ratio and large non-Gaussianity simultaneously. 

For the non-linearity parameters, we obtain 
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Since these expressions are rather complicated in general, we consider three limiting 
cases, which may generate large non-Gaussianity: Single-source case, equal mass 
case (?7i = r] 2 = rj), and large mass ratio case (771 3> 772). 



I. Single source case 

The setup of the multi-brid inflation model includes the inhomogeneous end of 
hybrid inflation discussed in section 13.2.41 as a limiting case where fluctuations 
of a single field, 0i or 02, only contribute to the curvature fluctuation and the 
angle a = 0. Here we regard 0i as the inflaton and 02 as another light scalar 
field, denoted as a in Section 13.2.41 in which the mass of the inflaton is much 
bigger than that of a light scalar field (i.e. 771 ^> 772). In this limit, we also 
assume that the inflaton does not contribute to the power spectrum, which is 
represented by the following condition: 

A 2 A 2 ( \ 

gffiljfa < S202J01 < 2 2 2 02,/J • ( 174 ) 

The last inequality applies only when 771^ < 1, which leads to the relation 
0i,//0i — 02, //02- In order to realize large non-Gaussianity, we further assume 
that A\ <C A?,, which is equivalent to consider the situation where #202/ *C 
g 2 4>lf — v 2 . Putting these conditions to the above general expressions with 
a = 0, we obtain the following formulae, 
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TNL ^ ( q/nL ) , (177) 



where the following equations have been used: 



^ , G{i) ~ -2^ • (178) 
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These results reproduce those obtained for the inhomogeneous end of hybrid 
inflation model discussed in Section 13.2.41 and again notice that 5nl//nl is 
suppressed by the slow-roll parameter rjx- 



II. Equal mass case = r] 2 = r]) 

In the equal mass limit (771 = 772 = 77), a-dependence disappears [93] because of 
the symmetry. Hence, we can set a = without loss of generality. In such a 
case, we have 
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Then, the tensor-to-scalar ratio and the spectral index are reduced to 
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ije vNk . The non-linearity parameters in this limit 
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After some algebra, we can find the simple relation between /nl and <?nl: 

10 , 50 2 
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As is the same with the single-source case, the ratio of the non-linearity param- 
eters Pnl//nl is suppressed by the slow- roll parameter i] in this well. 
We also obtain the following relation, 
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where the ratio R is given by 
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from which we can calculate the size of tnl relative to /^ L . 
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Figure 5: Contours of /nl (red lines) and tnl (blue lines) on the 7-51/(72 plane for the 
equal mass case. The blue lines are for / NL = 5 (thick solid), 10 (dashed), 50 (dotted) 
and 100 (thin solid). The red lines are for r NL = 36/25 x 5 2 (thick solid), 36/25 x 10 2 
(dashed), 36/25 x 50 2 (dotted) and 36/25 x 100 2 (thin solid). We set Vl = 7] 2 = r] = 0.01 
here. From this figure, we can find that the upper right region on this plane corresponds 
to the single-source case because the equality tnl = 36/25 x is satisfied, which holds 
in the single source case. 

In Fig. El we depict contours of /nl (red lines) and tnl (blue lines) on the 
7-51/52 plane. As seen from this figure, we can find that the upper right region 
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Figure 6: /nl^tnl diagrams for the equal mass case. In the left panel, the red solid line is 
for the case with | — 7 = 0.01 and g\jg 2 being varied. In the right panel, the red solid line 
is for the case with g\j g 2 = 10 2,5 and 7 being varied. In both the panels, the blue dashed 
lines indicate the equality tnl = 36/25 x /^ L . We set 7/1 = r] 2 = r\ = 0.01 and a = here. 



on this plane corresponds to the single source case because the equality tnl = 
36/25 x /nl, which holds in the single-source case, is approximated satisfied. 
In Fig. [6l we also show the /nl~7nl diagrams in this limiting case. In the left 
panel, the red solid line is for the case with f — 7 = 0.01. In the right panel, 
the red solid line is for the case with gi/g 2 = 10 2 ' 5 . In both the panels, the blue 
dashed lines show the equality tnl = 36/25 x /^ L . 

III. Large mass ratio case (771 ^> 772) with A 2 <C 1 



In order to study the case of large mass ratio, we have set A 2 = in the 
expressions for the non-linearity parameters Eqs. f ll7ip -( TT73l) . However, we 
would like to stress that, as already pointed out in Ref. [93j, our expressions 
apply for an arbitrary mass ratio as long as A 2 = 0. From the expression for 5N 
given in Eq. f)165p . we can easily find that the condition A 2 = corresponds 
to neglecting the linear contribution from ^7. In such a case, the non-linearity 
parameters are given by 
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From these equations, we find that, even if A\ is order of unity, 



4 1_ 

A\ 

large mass 



ratio 771 S> 772 can generate the large non-Gaussianity. For large /nl, we can 
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approximately obtain the following relations, 




(192) 



(193) 



(194) 



Here, it should be noticed that #nl is of the order of and can become 
relatively large. This is because both fields significantly contribute to the non- 
linearity parameters as well as the power spectrum despite their large mass 
ratio. 

It should also be noted that, in the single source case discussed earlier, we 
have assumed not only the large mass ratio but also the large ratio of the field 
values, 0i and 02, in order to get the large non-Gaussianity. In such 
the curvature perturbation was effectively generated only from a single source 
at all orders. 

On the other hand, here, we have assumed only A 2 = 0. From Eq. (11661) . 
we can find that this case is different from the single source case because the 
curvature perturbation can be generated from multi-source at the second (and 
also the third) order. Hence, the above expressions of non-linearity parameters 
are different from those in the pure single source case. 



B. Linear potential model 

Next let us briefly discuss the linear multi-brid model whose potential is given by 



MM 
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where 



[g\ (0i cos a + 02 sin a) 2 + g\ (-0i sin a + 2 cos a) 2 ] x 2 + t ( X 2 
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(196) 



For this potential, the slow-roll equation of motion is 




(197) 
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Hence the total e-folding number is evaluated as 

N = — (0i - 0ij) , (198) 
mi 

where 0ij is the value of the scalar field at the end of inflation. As in a usual hybrid 
inflation model, the inflation ends when the following relation is satisfied: 

2 2/ \ 2 2 / / \2 / \ 

v — 9i (0i,/ cos a + 02,/ sin a) + g 2 (— 0i,/ sin a + 2) / cos a) . (199) 

We again parameterize 0ij and 02,/ as in Eqs. (I162p and (11631) . The field values at 
the end of inflation can be described in terms of the values of 0i and 02 at some 
time during inflation as 

— 01 - —02 = —0i,/ - —02,/ • (200) 
rrii m 2 m i ^2 

Perturbing the above equation, we have the following relation, 

1 r / 1 x , ( 1 90!,/ 1 d<j> 2J \ 

501 502 = 77^ 77^ #7 

nil m 2 \ m l 0*7 HT'2 Cry / 

+ 1 f 1 92 ^f _ J_&foA + I f 1 ^1,/ _ 1 d 3 2 ,/ \ 3. 

2 \mi c?7 2 m 2 c?7 2 / 6 ymj cfy 3 m 2 d'j 3 J 

(201) 

From Eq. (11981) . 5N is given by 

XKT 1 A,A 1 ^I./ X 1 1 X 2 1 1 <9 3 01,/ x 3 . on0 N 

<5iV = 50! H-^^T-c «3"" ' ■ ( 202 ) 

mi mi (77 2 nil oj* b nil o^ 6 

Substituting Eq. (120 ip into Eq. (I198p . we have the expression of 5N up to the third 
order, 

5N = J— (-l 2 (7)50i + + (™ 2 50i - ^i50 2 ) 2 



F( 7 ) v > 2F( 1 f 

m 2 50i - mi50 2 ) 3 , (203) 



2F{ 1 y 

where 

A(7) = ^, *(7) = ^, (204) 
F( 7 ) = m 3 Ai (7) - miA 3 ( 7 ) , G( 7 ) = ^(7)^(7) - ^(7)^(7) , (205) 
and we have used 

A'l = -A u A" = -A 2 . (206) 
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From these equations, the non-linearity parameters can be calculated as 
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By using the above expressions, we find that #nl can be written with /] 
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Thus, (7nl m the linear potential model is at least of the order of f^, and it can be 
relatively large compared to /nl- 

The spectral index and the tensor to scalar ratio are given by 

n.-U-m+n^, r= ~ ^ (211) 

from which we find another simple relation among /nl,^s and r: 

(6/nl/5)2 = 1 - , r , . (212) 
r NL 8(l-n s ) 

Although the expressions for the non-linearity parameters in this model are not so 
complicated, but still it is a bit difficult to see its size explicitly. Thus we look at 
some limiting cases in order. 

I. Equal mass case (mi = m 2 = M) 

Let us first consider the equal mass case. The expression of /nl m this case is 
given by 
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Regarding tnl, from Eq. f)212p we obtain the following relation 



1 + iA /6 



**L = ( g/NL] , (214) 



where the ratio i? is given by 

fl = fi±iV. (215) 

We can also find the relation between /nl and #nl a s 

<7 2 cos 2 7 + g\ sin 2 7 sin 7 (cos a + sin a) — g 2 cos 7 (cos a — sin a) 2 

0NL = 2 : 7 ; — : : ; ; r / NL 

0102 g2 sin 7 (cos a + sin a) — g\ cos 7 (cos a — sin a) 

0(W02) x /nl, 01 > 02 

(216) 

O(02/0l) X /jsjl, #2 > 01 

In particular, g NL = 2f^ L for the case with g\ = g 2 . 

II. Large mass ratio case (mi ^> m 2 ) with A 2 <C 1 

Now we briefly discuss the large mass ratio case. To study this case, we again 
set A 2 = as we did for the counterpart in the quadratic potential case. Then 
the expressions for the non-linearity parameters become 

6 mi mi0 2 ,/ , 917 x 

r/NL - — ~— , [HQ 
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™ = IttJU /nl J • (218) 

m 201,f\ ,2 



^ NL = 2 1-^^ / 2 L . (219) 

Thus, the linear potential model generically predicts large #nl relative to /nl and in 
some limiting cases, we have ~ /nl as we nave shown above. 

3.4 Constrained multi-source model 

In this section, we discuss a class of "constrained multi-source" model. In fact, models of 
this category are particularly related to so-called loop contributions. Thus we start with 
a general discussion including loop terms. 
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3.4.1 Expressions including loop contributions 

In the discussions up to the previous section, we have neglected loop contributions in 
the expressions for the power spectrum and non-linearity parameters. However, if we 
include such loop terms, a new type of model can appear. The expression for P^ and the 
non-linearity parameters, including one loop contributions, are respectively given by 

P c (k) = [N a N a + N ab N ab V 5 \n(kL)] P 6 (k) , (220) 
6, N a N b N ab + N a b N b c N c a VsHk ml L) 

t/nl = — - — — - , (221) 

5 (N a N- + N bc N^VsHhL)f 

N a N b N ac N c b + N a b N b c N c d N d a V s \n(k m2 L) 

Tnl = 5 i 

(N a N- + N^VsH^L)) 3 

where L is the size of the box in which perturbations are defined and we have intro- 
duced the notations k m \ = min{/c;} and k m 2 = min{|/cj + kj\,ki}. In deriving the above 
equations, we have truncated the expansion of Eq. flTJ at the second order in the field 
fluctuations. Strictly speaking, higher order terms contribute to the one-loop corrections. 
The full expressions for the one-loop correction are given in appendix [B] In this section, 
we do not take into account such higher order effects since it turns out to be negligible in 
many cases. Since the current observations indicate that the primordial curvature pertur- 
bation is almost Gaussian, the power spectrum should not be dominated by the one-loop 
contributions, from which we have the following constraint: 

TV A/" * TV TV a b 

1 > -if^Vs In(kL) ~ ab 2 P C \n(kL). (223) 
N C N C v ; (N C N C ) 2 C V ; V ! 

Now let us consider the relation between /nl and tnl in the case where the one-loop 
contributions dominate in the non-linearity parameters. In such a case, /nl and tnl are 
respectively given by 

jj/ NL = P ab M ab Vc, Hk ml L), (224) 
r NL = M ab M ab Vc, Hk m2 L), (225) 

where 

N h N C N h 

Pab = Mab = a c \ . (226) 

N d N d) (N d N d ) 2 1 ^ 
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Using the Cauchy-Schwarz inequality, we find 
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/nl) = (P ab M ab ) 2 V 2 ( \n 2 (k ml L) 
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ir 2 

TNL 



In (k ml L) , b . 



ln 2 (/c m iL) 

< rNL ln(fc m2 L)ln(fcL)' (22?) 

where in the final inequality we have used Eq. (I223p . We have also assumed that k is any 
wavenumber that satisfies In(kL) ~ \n(k m iL) ~ ln(/c m 2-L) = 0(1). Then we obtain an 
approximate inequality given by 

6 



tnl > ( q/nl ) , (228) 

up to the logarithmic corrections. Hence, even in the case where the one-loop contribution 
dominates in the non-linearity parameters, "the local-type inequality" Eq. (fT4"j) should be 
satisfied under the condition where the contribution from the one-loop term should be 
subdominant in power spectrum. 

3.4.2 Ungaussiton model 

In some cases, the first terms in the numerators of Eqs. (I22ip and (1222 j) are negligible com- 
pared to the second order terms. In such a case, the non-linear parameters are dominated 
by the one loop contributions. This kind of model has been discussed in [34ti36| [95] and 
called "ungaussiton" in (36] and "quadratic model" in (95]. Let us consider the simplest 
case where the curvature perturbation is given by 

C = iV0 + ^N a Ja 2 + • • • , (229) 

where is the inflaton and a is the so-called "ungaussiton." In this model, the expression 
for the non-linearity parameters given by Eqs. and ffTU]) cannot be adopted in this case 
any more and /nl, tnl and #nl are given by 

5/NL = jji , (230) 

TNL = JTq , (231) 

54 2>NlN aa N aaaa V a \og{k ml L) + a Nl a N aaa V a log(fc ml L) 

25^ nl = m ■ (232) 
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Since and N aaa are strongly suppressed in this model, the size of #nl would be very 
small. Thus the trispectrum is dominated by tnl- From the above expressions, we find 
the relation between /nl and tnl as 

Tnl = CP^ftH ~ 10 3 /^ 3 , (233) 

where we adopt the normalization for the power spectrum ~ 2 x 10 -9 and C is a 
constant defined as C = (6/5) 4//3 (ln(fc ml L)/[ln(fc m2 )] 4 / 3 ). This equation clearly contradicts 
with Eq. ( |228p if /nl is sufficiently large. However, it can be easily checked that, in such a 
case, the power spectrum is also dominated by the loop contribution, which does not satisfy 
our assumptions to derive Eq. ( I228p . Also notice that the curvature perturbation induced 
from a starts from the second order. Hence the fluctuations from a are completely non- 
Gaussian, which cannot explain the observed fluctuations. Thus we need another source 
of fluctuations to account for the observed almost Gaussian fluctuations, which could be 
those from the inflaton. In this sense, this model requires multi-sources in nature. This 
motivates us to call this kind of models "constrained multi-source model." 

In fact, when the initial value (the field value during inflation) of the scalar field er* is 
less than its quantum fluctuation during inflation Jer* ~ H in f/ (2tt), models of this category 
can be easily realized such as in the axion, the curvaton, the modulated reheating and so 
on. 



3.5 Other models 

In this paper, we have discussed various models generating large local-type non-Gaussianity 
focusing on models which can be accommodated in the framework of inflationary universe 
(fluctuations originating to quantum fluctuations during inflation) with adiabatic mode. 
However, there are other possibilities of generating large non-Gaussianity in some other 
frameworks and/or models where a simple parametrization of the non-linearity parameters 
may not be adopted. Here we briefly mention such other models of the local type. 

When one considers the generation of dark matter and baryon asymmetry of the Uni- 
verse, isocurvature perturbations can be generated in some situations. Since the shapes of 
the bispectrum and the trispectrum are different from the counterparts of the adiabatic 
ones, a simple parametrization of /nljJnl and #nl cannot be easily compared with the 



ones defined for adiabatic fluctuation s^ 21 ! . In fact, isocurvature fluctuations have already 
been severely constrained by cosmological observations such as CMB, however, the non- 
linear effect can be large in some cases, in particular, where fluctuations are dominated 
by the second order term. In this case, the situation is quite similar to the ungaussiton 
model discussed in the previous section. The issue of non-Gaussianity in models with 



# 21 However, we can have some approximate relation between /nl and / N l ■ While the former is defined 
in Eq. and has been discussed in this paper, the latter is defined to parametrize non-Gaussianity in 
models with isocurvature mode |95| . We can estimate how large non-Gaussianity from isocurvature 
fluctuations can be compared to the counterpart in adiabatic mode by using the relation between these 
non-linearity parameters derived in |95j . 
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isocurvature fluctuations has been investigated by several authors. We refer the reader to 
Refs. [951 - 1100] for the detailed discussion on this issue. 

Another local-type model is so-called the Ekpyrotic scenario [10111103] . in which the 
scale-invariant curvature perturbation and large non-Gaussianity can be generated. The 
non-linearity parameters /nl and #nl in this models have been calculated and the relation 
between / NL and #nl is obtained as <?nl ~ /nl [103] ■ 

There are also several studies on the possibility of generating large non-Gaussianity 
during preheating phase [104H113] . In this case, the super-horizon scale fluctuations might 
be strongly coupled to the evolution of small scale fluctuations. In Ref. |112j . the authors 
have introduced a new parameter F^l, which describes the non-linearity of the curvature 
perturbations and claimed that during preheating era there seems to be a possibility of 
generating large non-Gaussianity whose form is quite different from the usual local-type 
non-Gaussianity. 

Another possibility of local-type non-Gaussianity we would like to mention is a model 
in which the e-folding number is quite sensitive to the field value. In such a case, the 
truncation at the leading order or at the next leading order of Eq. will no longer give 
a correct answer. We then need to go to the higher order calculations until the series 
converge or to resort to the full order calculations. In appendix [Cj we give a simple toy 
model belonging to this class and briefly discuss some interesting features of the bispectrum 
and trispectrum. 

4 Discussion and Summary 

In this paper, we made a classification of models generating large local-type non-Gaussianity 
by using some consistency relations between the non-linearity parameters /nl, Tnl and #nl- 
The first key relation is the ratio of tnl/ (6/nl/5) 2 , by which we classify local-type models 
into three categories: 

• single- source model (tnl/(6/nl/5) 2 = 1) 

• multi-source model (tnl/(6/nl/5) 2 > 1) 

• constrained multi-source model (tnl oc f$ L ) 

We have also shown that the "local-type inequality" 



holds true not only for the case where the tree level contribution dominates in the non- 
linearity parameters but also for the case where the loop contribution dominates there. 
This inequality has been shown under the condition that a loop contribution is subdomi- 
nant in the power spectrum, which is required by current observations. To our knowledge, 




(234) 
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since all models generating local-type large non-Gaussianity known today should satisfy 
the "local-type inequality," if future observations confirm that this inequality does not 
hold, local-type models would be practically ruled out as a mechanism of generating large 
non-Gaussian primordial fluctuations. 

On the other hand, if future observations find large /nl of local type and probe the 
relation between tnl and /nl with some accuracy, satisfying the local-type inequality, we 
can see what category of models would be favored. However, even if we can pick up the 
one of these categories, as we have discussed, there still remain many possibilities for each. 
Thus we need a further classification to pin down the model of large non-Gaussianity. 
For this purpose, we can make use of the relation between /nl and #nl- We have shown 
that models can be further divided into three types according to the relative size of #nl 
compared to that of /nl as follows: 

• Suppressed #nl type (<?nl ~ [suppression factor] x / NL ) 

• Linear # NL type (# NL ~ /nl) 

• Enhanced <?nl type (#nl ~ /nl w hh n > 1 or n = 2 for many models) 

Thus if we further probe the relation /nl and #nl in future observations, we may find that 
only a few models survive by using the above categorizations. Then we can figure out what 
type of models are favored as a mechanism of the generation of primordial fluctuations. 

We have also worked out the above mentioned relations for various concrete models 
in this paper. Although models can be categorized rather rigorously by using the ratio 
t nl/ (6/nl5) 2 , the relation between / N l and #nl can significantly differ depending on some 
model parameters, in particular, in multi-source models. For example, let us take the 
multi-brid inflation model with quadratic potential, which was discussed in Section 13.3.31 
This model predicts that <?nl ~ ^/nl with being r\ a slow-roll parameter, which is of the 
suppressed #nl type, for the equal mass case, while the relation becomes #nl ~ f^ L , which 
is of the enhanced #nl type, for the large mass ratio case. In other words, the relation 
between /nl and #nl should be carefully investigated to discriminate a model because a 
model can predict quite different relations depending on its model parameters. However, 
it also means that the relation would be useful to explore the parameters of a model. 

In this paper, we have focused on non-Gaussianity in various models and not discussed 
much the tensor-to-scalar ratio or gravitational waves. In most models where large non- 
Gaussian primordial fluctuations are generated, the tensor-to-scalar ratio is considered 
to be very small in general. However, there are a few examples which can give a rela- 
tively large tensor-to-scalar ratio as well as generating large non-Gaussianity. One of such 
examples is a mixed model of inflaton fluctuations with some other source discussed in Sec- 
tion [33HJ In this model, the fluctuations from the inflaton can also be responsible for the 
curvature perturbation, and in such a case, the tensor-to-scalar ratio can be large, which 
could be detectable in near future observations. Another example is multi-brid inflation 
model discussed in Section 13.3.31 If both the tensor-to-scalar ratio and non-Gaussianity 
are found to be large in the future, above mentioned models may become the target of 
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detailed studies. This illustrates that a comprehensive investigation using non-Gaussianity 
and some other information such as gravitational waves can give much insight into the 
pursuit of the generation mechanism of primordial fluctuations. 

If three non-linearity parameters /nl, tnl and #nl are well determined in future obser- 
vations, we may be able to pin down the model of large non-Gaussianity and pick up a 
right model of generating primordial fluctuations. The classification by using the consis- 
tency relation among the above three parameters would be very useful to pursue the origin 
of the structure of the Universe and give a deep understanding of the early Universe. 
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Appendix 



A Full expression for ( in multi-curvaton model 

In Section |3.3.2[ we have discussed the multi-curvaton model, but investigated only some 
limiting cases. Here we give the full expression for the curvature perturbation in the model 
up to the third order. In general, the curvature perturbation ( can be written as 

C = CaCa(l) + CfeCb(l) + CaaCa(l) + ^bbCb(l) + CabCa(l)Cb(l) 

+C a aa(a(l) + C&6&C&(1) + C aa b(a(l)(b(l) + C a bbCa{\) Cb(l) • (235) 

We list full expressions for the coefficients such as C a and so on in the following. 
For the first order, 

n faljfal + 3)(/b2 - 1) / OQR x 

Ca - /al-3/ bl + 3 ' (236) 

q _ fal( — fblfb2 + fbl + fwp — 3(/&i — l)/&2 (237) 
fal — 3/f,l + 3 

For the second order, 

Caa = 777 77 7^T7 [/al(/&2 - 1) fal ~ %fal ( 2 /&l + fbl + 3 /b2 ~ 8) 

4(Jol - O/61 + 

-3/ Q 2 ! (2/£ + 4/ 2 2 + 12/ 62 - 13) + 9/ ol (5/ M - 2/ 2 2 - 6f b2 + 2) + 27(/ 61 - 1))] , (238) 
C » = JTf It TW2 [2/ai/ 6 2 i(/ & 2 - 1) + fli (-4/ 6 3 1 (/ 62 - 1) - 2/ 2 (/ 3 + 2/ 2 2 - 11/ M + 8) 

4 Ual ~~ -J/61 + -J J 

+/6i (4/ b 3 2 + 8/ 2 2 - 15/ 62 + 3) + f b2 (-2/ 2 2 - 4/ fe2 + 3)) 

-3/«i(/ M - 1) (2/&(/ M - 1) + Vli{fm - 1) + fbi (4/& + 8/ 2 2 - 15/ M + 3) 

-2/ M (2/ b 2 2 + 4/ 62 - 3)) - 9(/ 61 - l) 2 / b2 (2/ 2 2 + 4/ b2 - 3)] , (239) 

Cab = 77 77 — 7777 [faiihi - 1) ((9 - 2/ 2 x ) /^ + (/ ol + 3)/ w (/ 2 X - / ai (/ 2 2 + 3/ 62 - 5) 

[Jal — Ojbl + O) 

-3 (/ b 2 2 + 3/ b2 + 1)) - 3/ al / 6 3 ! + (U + 3) 2 / fe2 (/ fe2 + 3))] . (240) 
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For the third order, 



Caaa = ~T7T7~f A ^0,3 [/£(/«» " X ) M " Vtl (Wfcl + 3/ 6 | + 9/ fe2 - 37) 

+/ 4 (-12/2 + 2/m ( g/ 2 + 18/&2 _ 7? ) + 6/ 4 + 26/ 3 _ 48/ 2 _ 216/fe2 + 343 ) 

+3/* a (24/3 + 2/ 2 + g/fe2 _ 35 ) + /w ( 12/ 2 + 36/fe2 _ 23 ) 

+3 (6/ 4 + 26/ b 3 2 - ll/ 2 2 - 105/ b2 + 77)) + 9f 2 al (6f b \ + 2f bl + f b \ (6f b \ + 18/ 62 - 55) 

+f bl (-15/ 2 2 - 45/ fe2 + 128) + 3 (6/ fe 4 2 + 26f b \ + 3/ 2 2 - 63/ b2 + 15)) 

+27/ al (-9/ 3 + g/ 2 - 3/ 61 (6/ 2 2 + 18/ M - 23) + 6/ 4 + 26/ 3 2 + 15/ 2 2 - 27/ 62 - 17) 

-243(/ 61 - 1) (-2/ 61 + f b \ + 3/ 62 - 2))] , (241) 



a 



bhb 



[12/^1 (/ 6 2 - 1) 



24(/ a i - 3/ 61 + 3) 3 

+Vtlf 2 blUb2 - 1) ("24/2 _ ( 3/ 2 + g /fe2 _ g 7 j + 6/ 2 + 18/fe2 _ g) 

+2f al (-12/&(/ w - 1) + 2/ 4 (6/ 6 3 2 + 12/ 2 2 - 95/ b2 + 77) 
+A 3 i (6/ b 5 2 + 20/ 4 - 86/ b 3 2 _ ig2/ 2 + 595/fe2 _ 343 ) 
-3/w (6/ b 5 2 + 20/ 4 - 33/3 _ g6/ 2 + 126/fo2 _ 33) 

+6/ 6 i/ 62 (3/ b 4 2 + 10/3 - 4/ 2 2 - 18/ M + 9) - f 2 b2 (6/3 + 20/ 2 2 + f b2 - 18)) 

+6/ a 2 1 (/ bl - 1) (24/* (/«, - 1) + 2/ 4 (3/3 + 6/ 2 2 - 32/ b2 + 23) 

+3/ 6 3 i (4/ b 3 2 + 8/ 2 2 - 37/ b2 + 25) + 3f b \ (6/ fc 5 2 + 20/ 4 - 29/3 - 78/ 2 2 + 108/ 62 - 27) 

-12/ 6 i/ 62 (3/ b 4 + 10/3 - 4/ 2 2 - 18/ M + 9) + 3/ 2 2 (6/3 + 20/ 2 2 + f b2 - 18)) 

+9/ a i(/ 6 i - l) 2 (l2/ 6 5 1 (/ 62 - 1) + 28/ 4 (/ fe2 - 1) + 6/3(/ b2 - 1) (2/ 2 2 + 6/ b2 - 3) 

+Qf bl Ub2 - 1) (2/ b 2 2 + 6/ fe2 - 3) + 12/ 61 (/ 62 - l)/ 62 (3/3 + 13/ 2 2 + 9/ fe2 - 9) 

~Qfb2 (6fb2 + 20/ b 2 2 + f b2 - 18)) + 54(/ 61 - l) 3 / 2 2 (6/3 + 20/ 2 2 + f b2 - 18)] , (242) 



Caab = -rr? - 1 . ^ 3 [/ai(/w - 1) (54&/& + 9/ al (8/* + 2/ ol - 15) / 4 

4(/al - 3/ 6 i ~\~ 3) 

+ (/ai + 3) 3 f b2 (f b2 + 3) {2f al - 2/ al (3/ 2 2 + 4/ fe2 - 5) - 3) 

-3/ fe 3 ! (4/ 4 ! - 2/ a 3 x (3/ 2 2 + 9/ 62 - 35) - 3f 2 al (6/ 2 2 + 18/ M - 43) - 45/ al - 27) 

-(/ai + 3)/ 2 (24/ 4 - 2f al (6/ 2 2 + 18/w - 41) + 3f 2 al (2/ 2 2 + 6/ 62 - 55) 

+9/ai (9/ b 2 2 + 27/ b2 - 25) + 27 (/ 2 2 + 3/ fe2 + 2)) + (f al + 3) 2 f bl (6/ 4 

-VI (3/ b 2 2 + 9/ b2 - 19) + f aX (6/ 4 + 26/3 - 16/ 2 2 - 120/ b2 + 49) 

+3/ai (6/ fe 4 + 26/3 + 23/ 2 2 - 3/ b2 - 25) + 9 (2/ 2 2 + 6/ b2 + l)))] , (243) 
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Cabb - 777 77 TTvT [fal(fb2 - 1) (18 (4:fai + fal ~ 3) f bl 

+ 3) 2 /w (4/ fe 2 2 + 12/ 62 - 3) - f al (l2/ 4 + 52/ 6 3 2 + 13/ 2 _ 105/&2 + 15 ) 
-3 (12/4 + 52/ 3 + 46/ 2 _ 6/b2 _ 3)) _ 3 f b \ {4f al - 2f al (3/ 2 2 + 9f b2 - 35) 
-3/ai (6/ 2 2 + 18/ 62 - 35) - 18) + (-24/ 4 + 2j* (6/ 2 2 + 18/ 62 - 77) 
+3/ a 2 i (8/ 6 2 2 + 24/ fe2 - 29) - 9/ al (8/ 2 2 + 24/ fe2 - 45) - 27 (4/ 2 2 + 12/ 62 - 5)) 
Hfai + S)/ 2 ! (6/ a 4 x - 2f al (3/ 2 2 + 9/ b2 - 28) + f al (6/ b 4 2 + 26/' - 38/ 2 2 - 186/ fe2 + 163) 
+3/ al (12/ 4 + 52/3 + 15/ 2 _ 99/fe2 _ + g ( 6/ 4 + 26/ 3 + 2g/ 2 + mfb2 _ 8 ^ 

+54/ ol /& + (/ al + 3) 3 / b2 (6/ 6 3 2 + 26/ 2 2 + 21/ M - 9))] . (244) 



B Full expressions including one-loop corrections 



In subsection I3.4[ we have considered the one-loop terms in the expression for the power 
spectrum and the non-linear parameters, neglecting the third order terms in the SN for- 
malism. Here, we show the full expressions for the one-loop correction terms, including 
higher order ones in the SN formalism [99] . 

The expression for the power spectrum including the full one-loop correction terms is 
given by 

P c {k) = [N a N a + N ab N ab V s ln(fcL) + N a N a \V s ln(fc max L)] P s {k) , (245) 

where fc max is the cut-off scale of the power spectrum and the last term in principle can 
be removed by shifting the homogeneous value of the scalar field [30] . 

The non-linearity parameters including the one-loop corrections are given as follows. 
For /nl, we obtain 

7/nl = [N a N a + N ab N ab V 5 ln(fcL) + N a N ab b V s \n(k max L 
5 

V N a N b N ab + ( N ab N bc N c a + 2N a N bc N abc ) V s \n{k ml L) 

1 



f ( N a N ab N bc c + -N a N b N abc c \VsHk max L) 



(246) 



where k m \ = min{/cj}. Please notice that a term with fc max again can be removed by shifting 
the homogeneous value of the scalar field. Here we have assumed that \n(kiL) ~ ln(A; 2 L) ~ 
ln(k 3 L) and all of these are just represented as ln(fcjL), which enables us to factor out the 
dependence on such a factor in the power spectrum. Other two non-linearity parameters 
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are: 



[N a N a + N ab N ab VsHk t L) + N a N ab b V 5 \n(k max L 



)] 



-3 



TNL 



x N a N ab N bc N c + 2N a N ab N cd N bcd V s \n(k ml L) 

+ (N ab N bc N cd N da + N a N bc N b d N acd + iV a iV fc iV acd iV 6 cd ) V s Hk m2 L) 
(+N a N ab N bc N cd d + NaN^N^N^] V s ln(fc max L)l , 



(247) 
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[N a N a + N ab N°»V s ln(kiL) + N a N ab b V s Hk m ^L 



)] 



-3 



25 



x iV a iV b iV c jV abc + 3 (N a N bc N b d N acd + N a N b N cd N abcd ) V 5 Hk ml L) 
+ ( l -N a N b N c N abcd d + lN a N b N ahc N c A VsHk^L)} , 



(248) 



where k m2 = min{|A;j + kj\, ki}. 

C Full order calculation 

In this appendix, we briefly discuss the higher order perturbation effects on the bispectrum 
and the trispectrum. So far, our analysis has been based on the perturbative expansion 
of the e-folding number in terms of the field fluctuations. In order to calculate the power 
spectrum, the bispectrum, and the trispectrum, it was sufficient to consider the lowest 
order contribution to those quantities (or the next lowest order if the lowest one vanishes 
like in the ungaussiton case). This perturbative expansion is extremely accurate and 
powerful in many models including all the models discussed in this paper since the higher 
order contributions are highly suppressed. 

Yet, it is still logically possible to consider a model in which the e-folding number is 
quite sensitive to the field value. In such a case, the truncation at the leading order or 
at the next leading order will no longer give a correct answer. We then need to go to the 
higher order calculations until the series converge or to resort to the full order calculations. 

To see what happens if the higher order terms are taken into account, we give a simple 
toy model that allows full order calculations at some level. Here, by at some level, it means 
that the final expressions of the bispectrum and the trispectrum involve multidimensional 
integrals that defeats analytic calculation. Nevertheless, we can derive some interesting 
consequences from these results, as we will see below. 

Let us assume that the e-folding number depends on a field a by 




(249) 
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where (• • • ) denotes the spatial average. The magnitude of a field perturbation is not 
necessarily smaller than a. This type of model may be realized in the context of modulated 
reheating scenario if the coupling constant for the inflaton decay is multiplied by a factor 
like exp (— Then, using the 5N formula, the curvature perturbation is given by 



C(x) = exp 



5a{x) 



exp 



Sa(x] 



1 . 



(250) 



The full order calculation was also done in |114j . where the authors calculate the power 
spectrum of the field which depends on the Gaussian field by the sine function rather than 
the exponential one. 

By assuming that 5a is Gaussian, we write the two-point function in Fourier space as 

= WP(ki)S(h + k 2 ) . (251) 

We further assume the scale invariant power spectrum, i.e. P(k) = Hi/ {2k 3 ). Then, using 
the standard formula, 



exp 



exp 



d 3 x b(x)Sa(x) 



d 3 X\d 3 x 2 



d 3 q 
(2^ 



P(q)e 



iq-(x 1 -x 2 )l.( z! 



b(x 1 )b(x 2 ] 



(252) 



we find that the connected part of the three-point function of £ is given by 



(CfcjCjfcjjCjfcg/ 



d 3 Xid 3 x 2 d 3 x^ e 



iki -xi+ik2-X2+ikz-X3 



x exp 



d 3 q 



p{a){e iq 12 + e iq ' X2S + e iq ' X3V 



(253) 



where x 



Xi — Xj. This is non- vanishing only when k\ + k 2 + ks = 0. Then after a 
suitable change of integration variables, we can show that the bispectrum of £ is given by 



ij 



■B f (fci,fc 2 ,fca 



d x\d x 2 e 



x exp 



ik\-x\+ik2-X2 



d 3 q 



P(g)(e^ a?1 +e^ 2 + e^ 12 ] 



The momentum integration can be done analytically |114j . 



exp 



d 3 q 
(2^ 



P(q)e^ 3 



-) 

x J 



2d 



where L is the size of the box and (5 = H%/(8n a . ). Then the bispectrum becomes 



S c (/ci, k 2 , h 



d 3 x±d 3 x 2 e 



ik\-x\+ik2-x; 



2/3 / T \ 2/3 



(254) 



(255) 



(256) 
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It is not easy to do the integration over the real space coordinates. However, we can 
extract the scale dependence of the bispectrum. Indeed, under the scaling transformation 
ki — > Xki with A being constant, the bispectrum scales as 

B c (Xk h Xk 2 , Xk 3 ) = X~ 6+ ^B c (k u k 2 , k 3 ) . (257) 

Or, equivalently, the bispectrum runs with the scale as 

« = -6 + 6£. (258) 
dlogk 

The derivative should be understood that the shape of the triangle is kept unchanged. 
Therefore, the larger /3 is, the larger the bispectrum is on small scales. On the other hand, 
it is difficult to extract the shape dependence from Eq. (12561) and we stop pursuing it 
further. 

We can do the same for the four-point function. The trispectrum is given by 



T c {kuk 2 MM) = J d?x l d i x 2 d?x 3 e iki-s 1+ ik 2 -x 2 +ik s -x 3 



x (lY(lY(lY(±Y(±Y(±Y. (259) 



\XiJ \X 2 J \X 3 J \Xi 2 J \x 23 J \x 31 J 
Therefore, the trispectrum scales as 

T c (Xh, Xk 2 , Xk 3 , Xh) = X- 9+12(3 T c (h, k 2 , k 3 , h) , (260) 
and therefore, the trispectrum runs as 



d log 
dlogk 



-9 + 12/3. (261) 



Due to the same reason as the case of the ungaussiton, we assume that the leading 
contribution to the power spectrum comes from inflaton fluctuations. Therefore, 

^ssp 5 ' (262) 

where e is again the slow-roll parameter. From f !258j) and (I26ip . we find that a ratio 

2 

(263) 



B <Y /( T c\ _ B l 



P V I VI) p&' 

does not run with the scale. But it will in general depend on the shape of the triangle and 
the quadrilateral. 
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